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ABSTRACT 

We present an observation with the Chandra X-ray Observatory of the unusual radio source G359.23-0.82 
("the Mouse"), along with updated radio timing data from the Parkes radio telescope on the coincident young 
pulsar J1747-2958. We find that G359.23-0.82 is a very luminous X-ray source (L x [0.5-8.0 keV] = 5 x 
10 34 ergs s" 1 for a distance of 5 kpc), whose morphology consists of a bright head coincident with PSR J 1747- 
2958, plus a 45"-long narrow tail whose power-law spectrum steepens with distance from the pulsar. We 
thus confirm that G359.23-0.82 is a bow-shock pulsar wind nebula powered by PSR J1747-2958; the nebular 
stand-off distance implies that the pulsar is moving with a Mach number of ~ 60, suggesting a space velocity 
« 600 km s" 1 through gas of density w 0.3 cm . We combine the theory of ion-dominated pulsar winds 
with hydrodynamic simulations of pulsar bow shocks to show that a bright elongated X-ray and radio feature 
extending 10" behind the pulsar represents the surface of the wind termination shock. The X-ray and radio 
"trails" seen in other pulsar bow shocks may similarly represent the surface of the termination shock, rather 
than particles in the postshock flow as is usually argued. The tail of the Mouse contains two components: a 
relatively broad region seen only at radio wavelengths, and a narrow region seen in both radio and X-rays. We 
propose that the former represents material flowing from the wind shock ahead of the pulsar's motion, while the 
latter corresponds to more weakly magnetized material streaming from the backward termination shock. This 
study represents the first consistent attempt to apply our understanding of "Crab-like" nebulae to the growing 
group of bow shocks around high-velocity pulsars. 



Subject headings: ISM: individual: (G359. 23-0.82) - 
stars: winds, outflows 



1. INTRODUCTION 

Many isolated pulsars are observed to generate relativistic 
winds. The consequent interaction with surrounding material 
can generate a variety of complex, luminous, evolving struc- 
tures, collectively referred to as pulsar wind nebulae (PWNe). 
For several decades the Crab Nebula was regarded as the 
archetypal PWN, but in the last few years the realization has 
grown that PWNe can fall into a variety of classes, depending 
on the propertie s and evo lutionary state of the pulsar and of 
its surroundings (Gaensler 2004). 

Arguably the most spectacular such sources are pulsar 
bow shocks, which are PWNe confined by ram pressure 
owing to the pulsar's highly supersonic motion through 
surrounding material. Less than ten years ago, just a handful 
of such sources were known, predominantly seen through 
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Ha emission produced where the interstellar medium (ISM) 
was shocked by the pulsar's motion (Cordesi ri996l) . How- 
ever, recent efforts at optical, radio and X-ray wavelengths 
have identified m any n e w such PWNe and their pulsars 
(e.g.. iFrail etalJ I1996L Ivan Kerkwiik & Kulkarnj 1200 It 
Olber t et alJl2001t iJones et alJl2002t iGaensler et aljl2 002b ) ; 
see lChatteriee & CordesH2002t) and lGaensler et alJ J2004) for 
recent reviews. These results have consequently motivated 
renewed theoretical efforts to model these systems, through 
which one can infer the propertie s of pulsar winds and 
of the surrounding medium (e.g., [ Bucciantini & Bandiera 
1 200 11 iBucciantinl l2002t Ivan der Swaluw et alJ l2003t 
Chatteriee & Cordes 2004). Pulsar bow shocks are a partic- 
ularly promising tool in this regard, because they correspond 
to pulsars more representative of the general population than 
young pulsars like the Crab, and are unbiased, in situ, tracers 
of the undisturbed ISM. 

X-ray emitting electrons are a powerful probe because their 
short synchrotron lifetimes allow us to trace the current be- 
havior of the central pulsar, in contrast to the integrated spin- 
down probed by radio-emitting particles. It is thus not sur- 
prising that significant advances in understanding "Crab-like" 
PWNe have been made since the launch of the Chandra X- 
ray Observatory. This mission's high spatial resolution in the 
X-ray band has allowed us to make detailed s tudies of mag- 
netic fields and particles in pulsar winds (e.g.. Ga ensler et alJ 
2002a: iPavlov et al1l2003l) . While Chandra has successfully 
also identified X-ra y emission fr om a few pulsar bow shocks 
(Kasp i et alJl2001t lOlbert et al] 1200 It IStappers et all 12003). 
these sources have generally been very faint, so that little 
quantitative analysis has been possible. 
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1.1. G359.23-0.82: The Mouse 

Here we consider the most recently established instance of a 
pulsar bow shock, G359.23-0.82, also known as "the Mouse". 
G359.23-0.82 wa s discovered in a radio surve y of the Galac- 
tic Center region ( Yusef-Zadeh & Bally 1987), one of several 
unusual sources identified at that time in this part of the sky. 

A radio image of the system is shown in Figure ^ Clearly 
a "Mouse" is indeed an appropriate description: a compact 
"snout", a bulbous "body" and a remarkable long, narrow, tail 
are all apparent. The tail fades into the background ~ 12' 
west of the peak position of radio flux. Further to the west 
can be seen part of the rim of the supernova remnant (SNR) 
G359. 1-0.5. The bright head and orientation of the tail sug- 
gested that the Mouse is powered by an energetic source pos- 
sibly ejected at high velocity from the supernova explosion 
which formed SNR G359. 1-0.5. However, H I absorption ob- 
servations showed that G359.23-0.82 is at a maximum dis- 
tance of ~ 5.5 kpc , while SNR G359 .1-0.5 is an unrelated 
background object ( Uchida et al. 1992). 

The Mouse remained largely an enigma until 
iPredehl & Kulkarnil (U9951) used the ROSAT PSPC to 
show that the bright eastern end of this object was a source 
of X-rays. These archival ROSAT data are shown as contours 
in Figure ^ Two bright sources are seen to the immedi- 
ate southeast of the Mouse, SLX 1744-299 (upper) and 
SLX 1744-300 (lower), b oth thought to be X - ray binaries 
near the Galactic Center JSkinner et alJ 119871 ISidoli etafl 
12001 . Much fainter emission can be seen coincident 
with the head of t h e M ouse. Although the detection of 
IPredehl & Kulkarnil i 19951) lacked the resolution and sen- 
sitivity to carry out any detailed calculations, they argued 
that the radio and X-ray properties of G359.23-0.82 made it 
likely that this source was a bow-shock PWN, powered by 
a pulsar of spin-down luminosity £~2 x 10 36 ergs s" 1 and 
space velocity V ~ 400 km s" 1 . Subsequent deeper X-ray 
observations confirmed the power-law spectrum of X-ray 
emission from G 359.23-0.82, supporting this hypothesis 
(Sidol iet alll999h . 

Motivated by these results, we recently carried out a search 
for a radio pulsar associated with G359.23-0.82, using the 64- 
m Parkes radio telescope. We successfully identified a 98-ms 
pulsar, PSR J1747-2958, with a characteristic age t = 25 kyr 
and a spin-down luminosity E = 2.5 x 10 36 ergs s" 1 , coinci- 
dent with the Mouse's head (Cam ilo et al]l2002h . The po- 
sition of the pulsar as determined through pulsar timing by 
Camil o et alJ J2002I) is denoted by the small ellipse coinci- 
dent with the bright radio and X-ray emission at the Mouse's 
eastern tip. The high spin-down luminosity of the pulsar, its 
small characteristic age, and its spatial coincidence with the 
head of the Mouse, all argue that the pulsar is associated with 
and is powering the nebular emission. 

The Mouse thus is no longer a mystery, but rather presents 
itself as a spectacular example of an energetic, high veloc- 
ity, pulsar interacting with the surrounding medium. We 
here present a Chandra observation of G359.23-0.82, as 
well as updated radio timing parameters which confirm that 
PSR J1747-2958 is physically associated with this system. 
The high X-ray luminosity of the Mouse enables the first de- 
tailed study of a pulsar bow shock at high energies, and lets 
us begin building a link between the hydrodynamic behavior 
of ISM bow shocks and the relativistic properties of pulsar 
winds. We present our observations and data analysis in 
and our results in In ^|we carry out a detailed discussion 



of the structure and energetics of the X-ray bow shock driven 
by PSR J1747-2958, and compare these data with hydrody- 
namic simulations. We constrain properties of the ambient 
medium, of the shock structure around the pulsar, and of the 
nebular magnetic field, and consider what this bright PWN 
implies for the interpretation of other pulsar bow shocks seen 
in X-rays. 

2. OBSERVATIONS AND ANALYSIS 
2.1. X-ray Observations 

G359.23-0.82 and PSR J1747-2958 were observed with 
Chandra on 2002 Oct 23/24, using the back-illuminated S3 
chip of the ACIS detector jBurke et alJll997l) . The detector 
was operated in standard timed exposure mode, for which the 
frame time is 3.0 seconds. Data were subjected to standard 
processing by the Chandra X-ray Center (version number 
6.9.2), and were then analyzed using CI AO v3 . . 1 and cal- 
ibration set CALDB v2 . 2 3. No periods of high background 
or flaring were identified in the observation; the total usable 
exposure time after standard processing was 36293 seconds. 

The ACIS detector suffers from charge transfer inefficiency 
(CTI), induced by radiation damage. However, the effects of 
CTI for the back-illuminated chips are relatively minor, and 
have not been corrected for here. 

Energies below 0.5 keV and above 8.0 keV are dominated 
by counts from particles and from diffuse X-ray background 
emission. Unless otherwise noted, all further discussion only 
considers events in the energy range 0.5-8.0 keV. 

2.2. Radio Timing Observations 

PSR J1747-2958 was initially detected on 2 002 Feb 1 us- 
ing th e Parkes radio telescope, as reported by Cami lo et alJ 
(2002). We have continued to monitor this source at Parkes, 
and summarize here the timing observations through 2003 
Oct 6. 

The pulsar is observed approximately once a month, for 
~ 3 hr on each day, at a central frequency of 1374 MHz us- 
ing the central beam of a 13-beam receiver system. Total- 
power signals from 96 frequency channels spanning a band 
of 288 MHz for each of two polarizations are sampled at 1- 
ms intervals, 1-bit digitized, and recorded to magnetic tape 
for off-line analysis. Time samples from different frequency 
channels are added after being appropriately delayed to cor- 
rect for dispersive interstellar propagation according to the 
dispersion measure of the pulsar (DM = 101.5 cm" 3 pc). The 
time series is then folded at the predicted topocentric period 
(P ss 98.8 ms) to form a pulse profile. Each of these profiles is 
cross-correlated with a high signal-to-noise ratio profile cre- 
ated by the addition of many observations; together with the 
precisely known starting time of the observation, we obtain a 
time-of-arrival (TOA), and uncertainty, for a fiducial point (in 
effect, the peak of emission) of the pulse profile. 

PSR J1747-2958 is intrinsically very faint at radio wave- 
lengths, and the observed flux density fluctuates somewhat, 
likely owing to interstellar scintillation. The net result is that 
a handful of observations do not yield adequate pulse profiles 
and these are removed from subsequent analysis. We have re- 
tained a total of 26 TOAs, obtained from 82.6 hr of observing 
time, that we use to obtain a timing solution. 

3. RESULTS 
3.1. X-ray Imaging 
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The resulting X-ray image is shown in Figure [2] Apart 
from the Mouse (to be discussed below), the following other 
sources are apparent: 

(a) In the south-east corner of the field, the very bright 
X-ray binary SLX 1744-299, whose position in our Chandra 
data is (J2000) RA 17 h 47 m 25?90, Decl. -30°00'02."0. The 
image of SLX 1744-299 shows broad wings, resulting 
from both the point spread function (PSF) of the telescope 
and from dust scattering. The linear feature seen running 
approximately east-west through SLX 1744-299 is a read -out 
streak, produced by photons from this source hitting the CCD 
while charges are being shuffled across the detector. The core 
of the image of SLX 1744-299 suffers from severe pile up, 
resulting from two or more nearly contemporaneous photons 
being detected as a single event. No signal at all is detected 
in the central few pixels, because the total photon energy per 
frame exceeds the on-board rejection threshold of 15 keV. 

(b) 21 other unresolve d sources, de t ected using the 
WAVDETECT algorithm jFreeman et all 120021) . The 
most prominent of these is at coordinates (J2000) RA 
17 h 47 m 13?59, Decl. -29°59'16"9. This source contains 
1895 ± 44 counts in the energy range 0.5-8.0 keV, and 
lies 0"35 from the position o f the USNO-A2.0 star 0600- 
28725346 (Mo net etal1 IT998). We assume that this X-ray 
source is associated with this star, and that the offset between 
the X-ray and optical coordinates is an estimate of the uncer- 
tainty in positions determined from these Chandra data. (The 
error associated with the X-ray centroiding is negligibly small 
by comparison.) The corresponding positional uncertainty in 
each coordinate is therefore w 0."25. 

X-ray emission from the Mouse itself is also clearly de- 
tected, containing ~ 13000 counts in the 0.5-8.0 keV band, 
after applying a small correction for background. In Fig- 
ure |5J a) we show an X-ray image of the Mouse alone. This 
image has not had an exposure correction applied to it; over 
the entire extent of the Mouse the exposure varies by less 
than 1%. This image shows that the source overall has an 
axisymmetric "head-tail" morphology, with the main axis 
aligned east-west. To within the uncertainties o f the radio 
timing position as listed by Cami lo et al.l (|2002), the loca- 
tion of PSR J1747-2958 is coincident with the peak of the 
X-ray emission seen from the Mouse. (In ^3.5l below. we will 
present an improved position, with an uncertainty reduced by 
a factor of ^50 in each coordinate.) 

The brightness profile along the symmetry axis is shown in 
Figure |4] This demonstrates that to the east of the brightest 
emission, the vertically-averaged count-rate decreases very 
rapidly, dropping by a factor of 100 in just 3."5. Figure[3a) 
shows that this sharp fall-off in brightness forms a clear 
parabolic arc around the eastern edge of the source. Along 
the symmetry axis of this source, we estimate the separation 
between the peak emission and this sharp leading edge by re- 
binning the data using 0"12 x 0"12-pixels, and in the result- 
ing image determining the distance east of the peak by which 
the X-ray surface brightness falls by e" 2 = 0.14. Using this 
criterion, we find the separation between peak and edge to be 
1"0±0."2. 

Further east of this sharp cut-off in brightness, the X-ray 
emission is much fainter, but is still significantly above the 
background out to an extent 7" - 8" east of the peak. Fig- 
ure |3l a) shows that this faint emission surrounds the eastern 



perimeter of the source; this component can also be seen at 
lower resolution in Figure |2] In future discussion we refer to 
this region of low surface brightness as the "halo". 

To the west of the peak, the source is considerably more 
elongated. Figure |4] suggests that there are three regimes to 
the brightness profile: out to 4"5 west of the peak, a relatively 
sharp fall-off (although not as fast as to the east) is seen, over 
which the mean surface brightness decreases by a factor of 
10. Consideration of Figure |3j a) shows that this corresponds 
to a discrete bright core surrounding the peak emission, with 
approximate dimensions 5" x 6". We refer to this region as 
the "head". 

In the interval between 5" and 10" west of the peak, the 
brightness falls off more slowly with position, fading by a 
factor of ~ 2 from east to west. Examination of Figure |3j a) 
shows this region to be coincident with an elongated region 
sitting west of the "head", with a well-defined boundary. As- 
suming this region to be an ellipse and that part of this ellipse 
lies underneath the "head", the dimensions of this region are 
approximately 12" x 5". We refer to this region in future dis- 
cussion as the "tongue". 

The western edge of the "tongue" is marked by another drop 
in brightness by a factor of 2-3. Beyond this, the mean count- 
rate falls off still more slowly, showing no sharp edge, but 
rather eventually blending into the background > 45" west of 
the peak. Figure|3ja) shows that this region corresponds to an 
even more elongated, even fainter region trailing out behind 
the "head" and the "tongue". We refer to this region as the 
"tail". The tail has a relatively uniform width in the north- 
south direction of ~ 12", as shown by an X-ray brightness 
profile across the tail, indicated by the solid line in Figure [5] 
The tail shows no significant broadening or narrowing at any 
position. 

3.2. Comparison with Radio Imaging Data 

In Figure |3jb) we show a high-resolution radio image of 
G359.23-0.82, made from archival 4.8 GHz observations 
with the Very Large Array (VLA). This image has the same 
coordinates as the X-ra y image in Figure |3ja) (similar data 
were first presented by Yusef-Zadeh & Ballvl ll98^: an even 
highe r resolution image is presented in Fig. 2 of lCamilo et alJ 
120021) . The radio image shares the clear axisymmetry and 
cometary morphology of the X-ray data. The "head" region 
which we have identified in X-rays is clearly seen also in the 
radio image, in both cases showing a sudden drop-off in emis- 
sion to the east of the peak, with a slight elongation towards 
the west. The radio image shows a possible counterpart of 
the X-ray "tongue", in that it also shows a distinct, elongated, 
bright feature immediately to the west of the "head". How- 
ever, this region is less elongated and somewhat broader than 
that seen in X-rays. In the radio, the "tail" region appears to 
have two components, as indicated by the two contour lev- 
els drawn in Figure |3jb), and by the dashed line in Figure [5] 
Close to the symmetry axis, the radio tail has a bright compo- 
nent which has almost an identical morphology to the X-ray 
tail. Far from the axis, the radio tail is fainter, and broadens 
rapidly with increasing distance from the head. This compo- 
nent does not appear to fade significantly along its extent. The 
full extent of the radio tail can be seen in Figure [l] where it 
can be traced a further 12' to the west before fading into the 
background, along most of this extent h aving the appearance 
of a narrow , collimated tube (see also Yusef-Zadeh & Fjaiivl 
119871 11989I) . No counterpart to the X-ray "halo" is seen in 
these radio data. 
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3.3. X-ray Spectroscopy 

Spectra for G359.23-0.82 were extracted in six regions. 
The first five regions are defined by the annular regions ly- 
ing between five smoothed contour levels, as shown in Fig- 
ure [6] Regions 1 and 2 represent inner and outer regions of 
the "head", respectively; region 3 encompasses the "tongue"; 
regions 4 and 5 respectively correspond to the inner and outer 
parts of the X-ray "tail". In all five cases, only data en- 
closed by these contours which is to the west of RA (J2000) 
17 h 47 m 15."94 are considered, so as to avoid contamination by 
the spectrum of emission from the halo. Region 6 (not shown 
in Fig. [6) corresponds to the halo, and consists of data falling 
in an annular region centered on the X-ray peak, and lying 
between radii of 4"8 and 14." 1. Only data in this annulus ly- 
ing east of RA (J2000) 17 h 47 m 15."79 were considered, so as 
to avoid contamination by the spectrum of the head and other 
bright regions. 

For each of the regions under consideration, we computed 
the appropriate response matrix (RMF) and effective area file 
(ARF) using the CIAO script acisspec, which weights the 
RMFs and ARFs for well-calibrated 32 x 32-pixel sub-regions 
of the CCD by the flux of the source at that position, and then 
combines these to produce a single RMF and ARF for the 
region of interest. Once spectra were extracted, they were 
rebinned so that each new bin contained at least 30 counts. 
Spectra were subsequently analyzed using XSPEC vl 1.2.0. 

As can be seen in Figure|2j the CCD background in this ob- 
servation is dominated by the PSF wings and dust-scattered 
emission from SLX 1744-299. The background thus shows 
a significant spatial gradient across the field, determined by 
the distance from SLX 1744-299. To provide a good estimate 
of the background at the position of G359.23-0.82, we thus 
extract background counts within the annular region shown 
in Figure|2] corresponding to data lying between 2'.6 and 3'.9 
from SLX 1744-299, but excluding regions enclosing USNO- 
A2.0 0600-28725346, the read-out streak from SLX 1744- 
299, and G359.23-0.82 itself. This region contains 8754 ± 94 
counts within the energy range 0.5-10.0 keV, at an approxi- 
mate surface brightness of 0.2 counts arcsec" 2 . The spectrum 
in this region was subtracted from those in the six extraction 
regions of G359.23-0.82 in all future discussion. 

The data in regions 2-5 are all fit well by power-law spectra 
modified by photoelectric absorption: individual fits to each 
of these four regions result in spectra with foreground hydro- 
gen column densities Nh ~ (2-3) x 10 22 cm" 2 and photon 
indices T ~ 2. When one fits region 1 to a power law, one ob- 
tains a comparable absorbing column (Nh w 2.8 x 10 22 cm" 2 ) 
but a somewhat harder spectrum (r w 1.6). However, the fit 
is not particularly good (x 2 jv = 222/ '172 = 1 .29), mainly be- 
cause of a hard excess seen in data at energies 7-10 keV. This 
excess is not seen in fits to regions 2-5. Although the spec- 
trum of the background increases at these energies, region 1 
is much smaller and contains many more counts than regions 
2-5. It is thus unlikely that we are seeing the effects of back- 
ground in region 1 but not in regions 2-5. 

Rather, it seems likely that the X-ray emission in region 1 is 
of sufficient surface brightness to be suffering from pile-up, in 
which multiple low energ y photons are misconstrued as a sin- 
gle higher-energy event ( Davis 2001 ). Pile-up not only affects 
the spectrum of a source, but also produces an effect known 
as "grade migration", in which the charge patterns produced 
by adjacent photons are combined and reinterpreted as a sin- 
gle pattern, potentially different from the pattern produced by 



either incident photon. Generally, grade migration increases 
the possibility that a photon (or combination of multiple pho- 
tons) will be mistakenly identified as a particle background 
event. Thus an additional test for pile-up is to compare the 
number of "level-2" events (i.e., events with grades appropri- 
ate for X-ray photons) with the number of "level-1" events 
(i.e., all events telemetered by the spacecraft, regardless of 
grade). Sources which have surface brightnesses well above 
the background but which are free of pile-up will generally 
show a consistent ratio of count-rates in level-2 to level-1 
events. However, sources suffering from pile-up will show a 
reduced level-2/level- 1 ratio, because of grade migration. We 
have computed this ratio for the Mouse, and find that in re- 
gions 2-5, the ratio of level-2 to level-1 events is 0.97 ± 0.01 
(where the error quoted is the standard error in the mean). 
However, at the center of region 1 this ratio markedly drops 
to 0.88 ±0.01. This thus appears to be an additional signa- 
ture of pile-up at the peak of the X-rays from G359.23-0.82, 
confirming the inference made from the spectrum of region 1 . 

We have tried to account for the effects of pile-up on th e 
spectrum by incorporating the pile-up model of Davis ( 2001 ), 
the free parameters in this model be ing the "grad e morphing 
parameter" and the PSF fraction (see lDavisl2001l for details). 
Including this extra component in the spectral fit successfully 
accounts for the hard excess seen in the spectrum of region 1, 
and provides a good fit to a power law with best-fit parameters 
2Vff = (3.0±0.3)x 10 22 cm" 2 andT = 2.0±0.3. This spectrum 
is noticeably softer than the fit inferred without accounting for 
pile-up, as expected. 

To sensibly compare the spectra from all of regions 1-5, we 
assume that there is minimal spatial variation in the integrated 
hydrogen column density across the < 1' extent of the Mouse. 
We consequently have carried out a joint fit to the spectra of 
these five regions, requiring all five fits to have the same value 
of Nh, but otherwise leaving Nh as a free parameter. The 
photon index of each region was allowed to vary freely. The 
corresponding spectral parameters are listed in Tabled the 
spectra and model are plotted in Figure We find that this 
model is a good fit to the data, resulting in an absorbing col- 
umn Nh = (2.7 ± 0.1) X 10 22 cm" 2 and photon indices in the 
range 1.8 < T < 2.5. Most notably, there is a clear trend of 
an increasingly softer spectrum as one moves away from the 
"head", the outer region of the "tail" being steeper than the 
"head" by a factor Ar = 0.7 ± 0. 1 at 90% confidence. 

We caution that pile-up in region 1 has two additional ef- 
fects on the data. First, the incident count-rate is significantly 
higher than that detected, simply because the piled-up event 
list contains multiple events incorrectly interpreted as single 
events. Using the spectral parameters for region 1 listed in 
Table ^ we infer that in the absence of pile-up, the count- 
rate for region 1 would be 0.28 counts s" 1 in the energy range 
0.5-8.0 keV, 40% higher than that observed. This is in good 
agreement with the ef fects on count-rate due to pile-up es- 
timated in Figure 21 of Townslev et alJ d2002l) . Second, eve n 
when we account for pile-up using the model of Davis ( 2001 ), 
the inferred flux is still an underestimate. This is because 
the spectral fit does not incorporate those events incorrectly 
tagged as non-standard as a result of grade migration. The 
flux listed for region 1 in Tabled is thus only a lower limit; 
without carrying out detailed modeling of the effects of grade 
migration, it is difficult to estimate by how much the flux has 
been underestimated. 

We did not include region 6 (the halo) in the joint fit 
shown in FigureQ as it contains substantially less counts than 
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the other regions, and possibly has a more complex spectral 
shape. In the second half of Table [l] we list the parameters 
resulting from spectral fits to region 6 using both non-thermal 
(power law) and thermal (Raymond-Smith) models. (Consid- 
erably more sophisticated thermal models are available, but 
were not considered warranted given the low count-rate.) The 
best-fit thermal and non-thermal models are both good fits 
to the data, but result in estimates of the absorbing column 
which are inconsistent with those inferred from regions 1-5. 
We have also fit these models to region 6 with the absorption 
fixed at Nh = 2.7 x 10 22 cm" 2 . However, for both power-law 
and Raymond-Smith spectra, the corresponding fits are not 
good matches to the data. 

3.4. Spatial Modeling 

The Chandra image of the Mouse reveals that most of the 
X-ray flux from this source is contained in the bright "head" 
region, a close-up of which is shown in the left panel of Fig- 
ure [8] Of particular interest for interpreting this emission is 
to characterize what fraction, if any, of the "head" is in an 
unresolved source. In principle, one can answer this question 
by developing a spatial model for this region, and then de- 
convolving it by the PSF to determine the underlying spatial 
structure. 

Unfortunately, t he pi le-up which is present in the "head" 
(see discussion in 33.3l above) not only alters the spectrum in 
this region, but also distorts the shape of the PSF. Specifically, 
the effective count-rate in the core of the PSF is reduced as a 
result of pile-up, while that in the wings is mostly unaffected. 
Because of this effect, an image of an unresolved source will 
appears broader than the standard Chandra PSF. Since from 
Table [2 we know the incident spectrum of the piled-up emis- 
sion, we can account for this by simulating the expected piled- 
up PSF and fitting this to our image. However, as shown be- 
low, it is likely that the "head" region consists of a compact, 
piled-up, source sitting on top of a more extended component 
suffering from minimal pile-up. In this case, proper spatial 
modeling of this emission requires us to simultaneously de- 
convolve the compact component by the piled-up PSF, and 
the diffuse component by a PSF without pile-up. The shape 
of each PSF further depends on the incident spectrum, which 
likely is different for the compact and extended components. 
The required analysis would be extremely challenging, and 
is beyond the scope of this paper. To properly determine the 
spatial structure of this region, we plan future observations 
with the High Resolution Camera (HRC) aboard Chandra, 
which has slightly better spatial resolution than ACIS, and 
which does not suffer from pile-up. (The HRC also provides 
sufficient time-resolution to simultaneously search for X-ray 
pulsations from PSRJ1747-2958.) 

In the meantime, we here provide an illustrative example 
of a potential spatial decomposition of emission in the "head" 
region. Here we fit our model directly to the image, not tak- 
ing into account the effects of the PSF. We are thus unable 
to determine the true spatial extent of the underlying model 
components. 

We have carried out our analysis using the SHERPA fitting 
program within CI AO. Specifically, we extracted a 5" x 5" 
image centered on the "head", with 0.25 ACIS pixels = 0."123 
sampling. This region is completely contained within region 1 
as defined in Figure [6] Using SHERPA, we then created a 
model consisting of two gaussians plus a level offset. The 
center, amplitude, FWHM, eccentricity and position angle of 
each gaussian are all left as free parameters, as is the ampli- 



tude of the offset. 

This model was fit to the data using the Levenberg- 
Marquardt optimization method within SHERPA. The origi- 
nal image and best-fitting model are shown in the two left- 
most panels of Figure [8] The model matches the data well, 
as is demonstrated by the image in the center right panel of 
Figure|8] which shows that the residuals are of low amplitude 
and display no systematic structure. In the rightmost panel 
of Figure [8] we show the best-fit model at four times higher 
resolution, where the two gaussian components can better be 
distinguished. The parameters for this model and their uncer- 
tainties are listed in Table|2] 

To determine the expected extent of an unresolved source in 
the absence of pile-up, we have simulated the expected PSF 
using the Chandra Ray Tracer (CHaRT) 11 , for a source at the 
position of the peak of X-ray emission from G359.23-0.82, 
and having the incident spectrum listed for region 1 in Tabled 
Because this PSF is not piled-up, it is not directly comparable 
to our data, but it serves as a likely lower limit to the extent of 
any unresolved source. Using SHERPA to fit a gaussian to this 
PSF, we find that an unresolved, pile-up free source should 
have a FWHM of 0"95 ± 0."01 (with uncertainty quoted at 
90% confidence). 

The good fit of the data to the simple model shown in Fig- 
ure [8] makes it clear that the spatial structure of the "head" 
consists of at least two components at the resolution of Chan- 
dra. In our fit, one of these components is clearly extended, 
but is almost circular. The more compact component has 
a slightly larger extent than that expected for an unresolved 
source in the absence of pile-up, and shows a significant de- 
viation from circular symmetry. The orientation of the major 
axis of this component does not coincide with the east-west 
symmetry axis of the overall nebula shown in Figure 0a). 
Because of the complicated ways in which pile-up affects the 
shape of the PSF, we are unable to conclude from this anal- 
ysis whether this compact companion is slightly extended, or 
represents an unresolved source. As discussed above, future 
observations with Chandra HRC can resolve this issue. 

We can estimate the fluxes in the two gaussian compo- 
nents as follows. The total count-rates (0.5-8.0 keV) of the 
first (compact) and second (extended) gaussians are 0.037 and 
0.137 counts s -1 , respectively. We assume that both regions 
have the same spectrum as determined for all of region 1, i.e., 
a power law with absorbing column Nh = 2.7 x 10 22 cm -2 and 
photon index Y = 1.8. Using the response matrix and effec- 
tive area determined for region 1, we estimate that the unab- 
solved flux densities (0.5-8.0 keV) needed to reproduce these 
components are f x = 1.5 x 10~ 12 ergs cm -2 s -1 for the com- 
pact gaussian, and f x = 5.6 x 10~ 12 ergs cm" 2 s" 1 for the more 
extended component. Both these estimates probably underes- 
timate the true flux in these regions, because of the effects of 
pile-u p on count-rate and on grade migration as discussed in 
^3.3l above. 

As a final note, we comment that the sub-pixel imaging 
technique used to impro ve the sp atial resolution of some 
ACIS observations (e.g., iLi et alJ 120031) cannot be applied 
here. That method uses events which are split over several 
pixels to better localize the incident photon. However, in our 
case many multi-pixel events in the "head" are rather due 
to piled-up multiple photons, rather than split-pixel events 
from single photons. Sub-pixel imaging applied to this source 
would thus almost certainly produce misleading results. 

11 See http : / /cxc ■ harvard ■ edu/ chart/) . 
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3.5. Radio Timing 

Having obtained TOAs for PSR J1747-2958 as described in 
32.21 we have used the TEMPO 12 timing software to derive the 
pulsar ephemeris. TEMPO transforms the topocentric TOAs to 
the solar system barycenter and minimizes in a least-squares 
sense the difference between observed TOAs and those com- 
puted according to a Taylor series expansion of pulsar rota- 
tional phase. This fitting procedure returns updated model 
parameters (pulsar spin parameters and celestial coordinates) 
along with uncertainties and their respective covariances, as 
well as the post-fit timing residuals. 

Often in pulsar timing, especially for observations with 
poor signal-to-noise ratio, we find that the TOA uncertainties 
are significantly underestimated. This is reflected in a large 
nominal x\ l v for the fit, and is the case for our data. A com- 
mon remedy employed in order to ultimately obtain realistic 
parameter uncertainties is to multiply the nominal TOA un- 
certainties by an error factor so as to ensure xt/ p = 1 • We 
determined this factor using a segment of the data that is short 
enough in span so that no other unmodeled effects are visible 
(cf. below), and hereafter use TOA uncertainties increased by 
a factor of 3.4. 

Fitting the TOAs for PSR J1747-2958 in a straightforward 
manner with a pulsar model consisting of the usual minimal 
set of parameters (RA, Decl., rotation frequency v = l/P, and 
frequency derivative) yields residuals that are not featureless. 
These residuals appear to be due to rotational "timing noise", 
a common occurrence in youthful pulsars such as PSR J 1 747- 
2958. In the presence of timing noise, the parameters as esti- 
mated in this manner are biased, and we f ollow instead an al- 
ternat ive prescription for such cases (e.g., Arzoumani an et all 
1994). First, we fit for RA, Decl., pulse frequency, and as 
many of its derivatives as are required to absorb the unmod- 
eled noise and "whiten" the residuals. In our case this requires 
two derivatives. The resulting celestial coordinates and re- 
spective uncertainties represent our best unbiased estimates of 
these parameters, and are presented in Table[3] We then fix the 
pulsar position at the values thus obtained, and perform one fit 
for frequency and its two first derivatives. The resulting val- 
ues and uncertainties, given in Table[3] minimize the rms tim- 
ing residuals for this pulsar over the time span of the observa- 
tions. The values of v and v thus obtained are biased slightly 
when compared to their "deterministic" values (obtained by 
fixing v at zero) but the resulting ephemeris has better pre- 
dictive value near the present epoch. The value for the fre- 
quency second-derivative is unlikely to be deterministic, but 
rather gives information about the magnitude of timing noise: 
v measured for PSR J1747-2958 is consistent with a level of 
timing noise comparable to that experien ced by pulsars hav- 
ing si milarly large v (see compilation by Arzoumania n et all 
1994). 

Having followed the above procedure, we are confident that 
the celestial coordinates are free from significant systematic 
errors. The position of the pulsar thus o btained is consistent 
with that reported by ICamilo et alJ ( 120021) . based on a subset 
of the data we use here, but has much higher precision. Note 
that the uncertainty in Decl. in Table [3]is seven times larger 
than in RA, due to the low ecliptic latitude of the pulsar. 

Using the covariance matrix of the fit and the formal stan- 
dard errors returned by TEMPO for individual parameters, we 
can obtain for a given confidence level the joint confidence 

12 See |http : / /pulsar .princeton . edu/tempo| . 



region for RA and Decl. . The ellipse plotted in the first and 
fourth panels of Figure[8]displays this region for a confidence 
level of 99.73% (i.e., 3 a). Once one factors in the uncertainty 
of « 0."25 in each coordinate betw een t he radio and X-ray ref- 
erence frames (see discussion in 33. 1> . there is clearly good 
agreement between our updated position for the radio pulsar 
and the peak of the most compact region of X-ray emission. 

4. DISCUSSION 

On t he basis of a faint X-ray detection using the ROSAT 
PSPC. lPredehl & Kulkariill ill 9951) argued that the Mouse was 
a pulsar-powered bow shock. From this interpretation one 
could make two clear predictions: that higher resolution 
X-ray imaging would show a non-thermal cometary nebula 
(most likely of smaller extent than that seen at radio wave- 
lengths), and that an energetic radio pulsar would be found 
coincident with the head of the radio/X- ray nebula. 

The detection of PSR J1747-2958 bv lCamilo et ail (l20()l 
identified the likely central engine. Our Chandra image, com- 
bined with our updated radio timing data, confirms both that 
G359.23-0.82 has the expected X-ray morphology, and that 
PSR J1747-2958 is almost perfectly positionally coincident 
with the Mouse's bright head. We therefore conclude that 
G359.23-0.82 is indeed a PWN associated with PSR J1747- 
2958, the pulsar's supersonic motion from west to east gener- 
ating a bow-shock morphology. 

Other than PSR J1747-2958 and G359.23-0.82 we identify 
five other instances in which X-ray emission from a pulsar 
bow shock has been confirmed, as listed in Table|4] However, 
in the observation presented here, we detect at least an order 
of magnitude more X-ray photons from this source than seen 
so far from most of these other sources. 13 This system thus 
presents our best opportunity yet to study the X-ray emission 
produced by the interaction between a supersonic pulsar and 
its surroundings. 

Before considering the detailed structure of this source, we 
make one comment on its X-ray spectru m compared to o ther 
young pulsars and their PWNe. Recently, Gotthelf (2003) has 
argued that for many such systems, the photon index of the 
power-law X-ray spectrum seen for both pulsar and PWN are 
correlated with the pulsar's spin-down luminosity. Specifi- 
cally, pulsars of progressively lower values of E appear to 
have progressively flatte r photon indices. The least ener- 
getic pulsar presented by Gott heB (120031) is PSR J181 1-1925, 
which has a spin-d own luminosity E = 6.1 x 10 36 ergs s -1 
iGavriil et al J 120 04). with photon indices fo r the pulsar an d 
PWN of T = 0.6 and T = 1.3, respectively jGotfhelfl t2003). 
Clear ly the Mo use does not follow the trend described by 
iGotthelfl 12003): despite PSR J1747-2958 having a value of 
E a factor of > 2 lower than for PSR Jl 81 1-1925, the pho- 
ton indices listed in Tableware much steeper for this system 
than for PSR J1 81 1-1925. Th is is consistent with the con- 
clusion made by Gotthelf (2004), namely that the relationship 
described by Gotthelf (2003) does not extend to pulsars with 
E < 3.5 x 10 36 ergs s -1 , or to PWNe with bow-shock, rather 
than "Crab-like" morphologies. 

4.1. Distance to the Mouse 

13 The exception is PSR B1951+ 32, for which ~60 000 counts were de- 
tected in the Chandra observation of Moon et al. 12004). However, this sys- 
tem reflects a complicated interaction with the associated SNR CTB 80, and 
so is difficult to interpret. 
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The distance to G359.23-0.82 and to PSR J1747-2958 is 
something we reconsider here. The lack of H I absorption 
from G359.23-0.82 seen against the 3-kpc spir al arm argues 
that t he distance to this source is < 5.5 kpc (lUchida et alJ 
1992). The disp ersion measure tow ards the pulsar is DM = 
101.5 cm" 3 pc (Camiloetal. 2002), which implies a dis- 
tance of 2 kpc when usi n g the Galactic electron density 
model o f ICorde s & Lazio| (|2002). (The earlier model of 
Tavlor & Cordes 1993 yields a distance re 2.5 kpc.) 

Using the X-ray spectrum of G359.23-0.82, we have here 
measured an absorbing column to this system of Nh ~2.7x 
10 22 cm -2 (see Table E). Assuming that G359.23-0.82 and 
PSR J1747-2958 are associated (which seems virtually cer- 
tain given their almost exact spatial coincidence), this implies 
a ratio of neutral hydrogen atoms to free electrons along the 
line of sight Nh/DM = 85. This is a value higher than that 
seen for all other X-ray detect ed pulsars, for whic h typically 
we observe 7%/DM re 5 - 10 l lDib & Kaspill2004l) . We pro- 
pose that this high ratio of neutral atoms to electrons can 
be explained if the Mouse lies behind significant amounts of 
dense molecular material. Indeed the only other pulsars with 
comparable ratios o f ' N H fDM are PSR B 1853+01, for which 
Nh /DM re 30-70 dWolszczan et al.lll99ll iRho et al.lU99H 
Harm s et all 19971) and PSR B 1757-24, which h as Afo/DM re 
40 (lManchesteretal.lll991t iKaspi et all 1200 ll) . It has long 
been established that PSR B 1853+0 1 and its SNR ar e embed- 
ded in a molecular cloud complex ( Seta et aLl l2004l and ref- 
erences therein), while PSR B 1757-24 lies behind the well- 
know n molecular ring at Galactocentric radius ^3-5 kpc 
(e.g.. lClemens et al.ll9 88). The Mouse lies just six degrees on 
the sky from PSR B 1757-24, and is similarly aligned in pro- 
jection with the molecular ring. However, the molecular ring 
is more distant from the S un than the distance to the Mouse of 
2 kpc adopted by Camilo et al.l (120021 ). To reconcile the high 
value of Nh seen here, we suggest that the Mouse is at a larger 
distance than 2 kpc, lying in or behind the molecular ring. 

To quantify this suggestion, we have integrated along the 
line of sight the equivalent hydrogen column density due to 
molecular clouds, using the radial profile of molecular sur - 
face density in the Galaxy given in Figure 1 of D amd d!993h . 
We assume that the m olecular layer has a FWHM of 120 pc 
jBronfman et alJl988l) 14 : because the Mouse is at b = -0.°8, at 
successively greater distances from Earth its sight line inter- 
sects molecular material increasingly further from the Galac- 
tic plane, and thus passes through material of relatively lower 
density. 

When we perform this calculation, we find that if the dis- 
tance to the source is 2 kpc as suggested by the pulsar's 
DM, the equivalent column due to molecular gas is only 
Nh ~ 4 x 10 21 cm" 2 . Unless the pulsar fortuitously lies behind 
a local molecular cloud, it is difficult to reconcile a distance 
of 2 kpc with our observed value of Nh- 

However, if the distance to the source is 4-5 kpc, signifi- 
cant amounts of dense material from the molecular ring lie in 
the foreground to this source, and the estimated column due 
to molecular gas increases to Nh re (1.0— 1.3) x 10 22 cm" 2 . 
Assuming that approximately half of neutral gas by mass is 
in atomic rather than molecular form (Dame 1993), we there- 
fore propose that the distance to the Mouse is approximately 
5 kpc. This is consistent with both the observed total column 
Nh = 2.1 x 10 22 cm" 2 , as well as with the upper limit to the 

14 We have scaled the results of Bronfman et al. 1 1988) to a Galactocentric 
radius of 8.5 kpc. 



distance inferred from H I absorption bv lUchida et al.l ( fl992l) . 
While this distance is in disagreement with that implied by the 
pulsar's DM, it is reasonable to expect larg e uncertainties in 
the Galactic electron density model of Corde s"& Laziol (120021) 
in this direction, as there are few known pulsars towards the 
Galactic Center to calibrate the distribution. At a distance of 
5 kpc, the pulsar's DM implies a mean free electron density 
along the line of sight n e = 0.0 20 cm" 3 , a not atypic al value 
towards the inner Galaxy (e.g., Johnst on et alJl20Ql|) . In fu- 
ture discussion, we therefore assume a distance d = 5ds kpc, 
arguing from the above discussion that d$ re 1 . 

We note that even at this increased distance, an association 
of the Mouse with SNR G359. 1-0.5, from which it appears 
to be emerging in projection, is still unlikely. Hi absorp- 
tion against G359. 1-0.5 clearly shows it to be at a distance 
of 8-10 kpc (inchida et al1 l!992l) . still significantly more dis- 
tant than the Mouse. This is further supported by X-ray ob- 
servations of G359. 1-0.5, which i mply an absorbing c olumn 
Nh ~ 6 x 10 22 cm" 2 to this source (Bamba et al. 2000), more 
than double that seen here towards G359.23-0.82. 

The revised distance to G359.23-0.82 immediately sets 
the size and luminosity of the system. The full radio ex- 
tent of the Mouse seen in Figure is an incredible 17^ pc, 
while the X-ray tail has length ~ I.W5 pc. The total unab- 
solved isotropic X-ray luminosity in the range 0.5-8.0 keV is 
Lx = 5dj x 10 34 ergs s" 1 . This is 2d\% of the pulsar's total 
spin-down luminosity, a conversion efficien cy exceeded only 
by a few other pulsars dPossenti et a l. 2002), and much more 
than the low values of Lx IE seen for other X-ray pulsar bow 
shocks, as listed in Table 0] Even if one adopts a nearer dis- 
tance d$ re 0.4 as argued by Cam ilo et al.l (|2002), the X-ray 
efficiency of the Mouse is still significantly above these other 
systems. 

It has previously been noted by several authors that bow- 
shock PWNe are particularly inefficient at converting their 
spin-down luminosity into X-ray emission; possible expla- 
nations i nclude minimal synchrotron cooling in the emitting 
regions (ICheyalie rl 120001) , parti cle acceleration in only a lo- 
calized area (Kaspi et al. 200l|, or weakening of the termi- 
nation shock because of the signifi cant mass-load ing present 
in ram-pressure confined winds (Lvutikov 2003). However, 
the Mouse demonstrates, as is also the case for "Crab-like" 
PWNe, that the X-ray efficiencies of bow shock PWNe can 
span a wide range, even for pulsars of similar age and 
E. Clearly environment, evolutionary history, magnetic field 
strength and possibly orientation of the system with respect to 
the line of sight all play an important role. 

4.2. Compact Emission Near the Pulsar 

We have shown in S13.4I and Figure [8] that the bright X-ray 
emission in region 1 can be well-modeled by two gaussian 
components, of FWHMs 1"1 and 2 "4 respectively. While the 
former is slightly larger t han that expected for an unresolved 
source, we have noted in 33.4l that this apparent extension may 
be entirely due to pile-up in this source. While future obser- 
vations can properly address this issue, for now we assume 
that the more compact gaussian represents an unresolved but 
piled-up source. 

In this case, given this source's location near the apex of 
the PWN, plus the good match of this location to the radio 
timing position of PSR J1747-2958, the most likely explana- 
tion is that this compact component represents X-ray emis- 
sion from the pulsar itself. The flux densit y wh ich we in- 
ferred for the first gaussian component in 33.41 implies an 
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isotropic X-ray luminosity (0.5-8.0 keV) for the pulsar of 
L x > 4.5fi?| x 10 33 ergs s" 1 . 15 This represents > 10% of the 
total X-ray emission produced by this system, a value typical 
of other young pulsars and their PWNe. 

In Table ^ we have only listed the fit to this spectrum 
from a power-law model. Fits of poorer quality can be made 
to blackbody models, but the inferred fit parameters (Nh ~ 
1 .7 x 10 22 cm -2 , kT w 1 keV and R BB w 200 m, where R BB is 
the equivalent radius of an isotropic spherical emitter) are not 
consistent with the column density of this source seen for re- 
gions 2-5, or with the expected properties of thermal emission 
from the surface of young neutron stars. We thus think it most 
likely that the emission from this source is non-thermal, and 
that it represents emission from the pulsar magnetosphere, 
which should be strongly modulated at the neutron star ro- 
tation period of 98 ms. Future X-ray observations of higher 
time resolution can confirm this prediction. 

Figure|8]demonstrates that in addition to the compact com- 
ponent which we have associated with PSR J1747-2958, the 
"head" also can be decomposed into a second extended gaus- 
sian, close to circular and 2.4" in extent. This component is 
by far the single most luminous discrete X-ray feature seen 
within the Mouse, with an X-ray luminosity (0.5-8.0 keV) 
of L x = Udj x 10 34 ergs s" 1 , more than 30% of the total 
X-ray flux from this system. As we will argue in 34.51 be- 
low, this region is contained entirely within the termination 
shock of the pulsar wind, in the region where the wind is 
cold and generally not generating any observable emission. 
However, for PWNe associated with both the Crab pulsar 
and with PSR B 1509-58, observations with Chandra have 
revealed compact X-ray knots at comparable distances from 
the pulsar, located within the wind termination shock and 
with spect ra harder than for the rest of the PWN, j ust as 
seen here ( Weisskoo f et all2000UGaensler et al]l2002al) . The 
origin of these knots is not known, but they ma y represent 
quasi-stationary shoc ks in the free-flowing wind (iLoull 19981 
Gaensler et al. 2002a). Although lack of spatial resolution 
prevents us from saying anything definitive here, we speculate 
that the bright extended region of X-rays seen here immedi- 
ately adjacent to PSR J1747-2958 may represent such knot 
structures produced close to the pulsar. Certainly G359.23- 
0.82 would then be unusual in having such structures repre- 
sent such a large fraction of the total X-ray luminosity. Sig- 
nificant Doppler boosting in this highly relativistic flow may 
be a possible explanation for this. 

4.3. Bow Shock Structure and Contact Discontinuity 

To interpret the other structures we see in the emission from 
G359.23-0.82, we have carried out a hydrodynamic simula- 
tion to which we can compare our data. In subsequent discus- 
sion we will show that the input parameters chosen for this 
simulation are a reasonable match to those likely to be appli- 
cable to the Mouse. We note that several previous simulations 
of pulsar bow shocks exist in the literature; those carried out 
by Bucciantini (2002) and van der Sw aluw et al.l (120031) are 
most rele vant to the data con sidered here. However, the sim- 
ulation of iBucciantinil (1200 2) does not i ncorporate regions far 
behind the apex of the bow shock, while lvan der Swaluw et al.l 
(2003) did simulate regions significantly far downstream, but 
only considered pulsars with low Mach numbers through the 
ISM. Our new simulation incorporates shocked material far 

15 As emphasized earlier, this value is only a lower limit because of the 
effects of pile-up. 



from the apex and considers high Mach numbers, both of 
which are likely to be relevant for understanding the Mouse 
(see further discussion below). 

We proceed in t h e sam e manner as described by 
Ivan der Swaluw et al.l (|2003) in their bow shock simula- 
tions. We simulate a pulsar moving at a velocity of 
600 km s" 1 through a uniform medium of mass density p = 
7 x 10~ 25 g cm" 3 . For cosmic abundances we can write 
p = l.31noniH, where nin is the mass of a hydrogen atom 
and «o is the number density of the ambient ISM; the density 
adopted thus corresponds to a number density «o = 0-3 cm" 3 . 
We used the Versatile Advection Code (VAC) 16 ilTothll 19961) 
to solve the equations of gas dynamics with axial symmetry, 
using a cylindrical coordinate system in the rest frame of the 
pulsar. A non-uniform grid was adopted, centered around the 
pulsar so that the pulsar wind could be resolved. The to- 
tal number of grid cells was 150 x 150. We used a shock- 
capturing, Total-Variation-Diminishing Lax-Friedrich scheme 
to solve the equations of gas dynamics over the total grid. This 
scheme yields a thickness for both shocks and contact discon- 
tinuities of typically ~ 4 grid cells. 

Apart from the differences in the Mach Number, Ai, and 
the wind luminosity, E, of the pulsar, th e simulation is i nitial- 
ized in the same way as described bv Ivan der Swaluw et alJ 
( 120031) : the pulsar wind is simulated by depositing energy at 
a rate E = 2.5 x 10 36 ergs s" 1 (as is observed for PSR J1747- 
2958) and mass at a rateM = 5.56 x 10 17 g s" 1 continuously in 
a few grid cells concentrated around the position of the pul- 
sar. The terminal velocity of the pulsar wind is determined 
from these two parameters, i.e., = (2E/M) X I 2 = 0.1c, and 
has a value much larger than all the other velocities of interest 
in the simulation. The pulsar wind velocity converges toward 
the terminal velocity; the associated Mach number has a max- 
imum value of ~ 20. 

This initialization of the pulsar wind results in a roughly 
spherically symmetric pulsar wind distribution before it is ter- 
minated by the surrounding medium. The current version of 
the VAC code does not include relativistic hydrodynamics, 
therefore for simplicity and for accuracy in the non-relativistic 
part of the flow, we adopt an adiabatic index for both rela- 
tivistic and non-relativistic fluid 7 = 5/3; we defer a full rel- 
ativistic simulation to a future study. 17 It is expected that a 
simulation which would treat the pulsar wind material as rela- 
tivistic and the ISM material as non-relativistic would slightly 
change the stand-off distance of the pulsar wind, but would 
not qualitatively change the results of the current sim ulati on. 
We adopt a Mach number M. = 60, which we show in 34.4l be- 
low likely describes the situation for the Mouse. This is much 
lar ger than the value M. = 7/V 5 « 3 . 1 , which was adopted 
bv lvan der Swaluw et al.l (120031) for the case of a pulsar prop- 
agating through a supernova remnant rather than through the 
ambient ISM. 

We perform the simulation using the abovementioned pa- 
rameters until the system is steady. We then multiply the 
length scales in the final output by a factor 10 -1 / 2 , because the 
terminal velocity used in the simulation (Voo = 0.1c) yields a 
stand-off distance 10 1 ' 2 larger then for a more physical ter- 

16 See |httpT/ /www .phys . uu . nl/ ~toth7| . 

17 We note that Bucciantini 1 2002) has carried out a bow shock simula- 
tion in which he distinguished between relativistic and non-relativistic ma- 
terial by adopting 7 = 4/3 and 7 = 5/3 for these two components, respec- 
tively. However this approach produced similar results to those obtained by 
van der Swaluw et al. 1 2003) using a uniform index 7 = 5/3. 
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minal velocity of Voo = c (see Equation (2) below). 18 The re- 
sulting bow shock morphology is depicted in Figure[9] which 
shows a logarithmic gray-scale representation of the density 
distribution. The scale of features in the simulation should be 
directly comparable to the data. 

This simulation clearly reveals the multiple zones and in- 
terfaces seen in previous simulations both of pulsar bow 
shocks i lBucciantinil 120021: Ivan der S waluw et al.| |200 3|) and 
of other supersonic systems (e.g., iMac Low et al.l 1199 It 
IComeron & Kaperl fT998t iLinde et alJ 119981) . Moving out- 
wards from the pulsar, these regions are as follows: 

A. Pulsar Wind Cavity: Immediately surrounding the pulsar 

is a region in which the relativistic wind flows freely 
outwards. Particles in this wind are assumed to have 
zero pitch angle and to not produce significant emis- 
sion. This unshocked wind zone is yet to be observed in 
any bow shock, but in X-ray images of the Crab Nebula 
and other Crab-like PWNe, this region is clearly visible 
as an underl uminous zone immediately s urrounding the 
pulsar (e.g.. IWeisskopf et alJl2000l iLu et alJl2002l) . At 
the point where the energy density of the pulsar wind 
is balanced by external pressure, a termination shock 
(TS) is formed where particles are thermalized and ac- 
celerated. In the Crab Nebula and other related sources, 
this interface is seen as a bright ring or arc surrounding 
the underluminous zone. Figure|9]shows that for a bow 
shock, the TS is highly elongated, having a significantly 
larger separation from the pulsar at the rear than in the 
direction of motion. 

B. Shocked Pulsar Wind Material: Beyond the TS, parti- 

cles gyrate in the ambient magnetic field and gener- 
ate synchrotron emission seen in radio and in X-rays. 
There are two distinct regions of emission in this zone. 
The flow near the head of the bow shock advects the 
synchrotron emitting particles back along the direction 
of motion of the pulsar, yielding a broad cometary mor- 
phology marked as region Bl in Figure [9] Directly be- 
hind the pulsar, material shocked at the TS flows in a 
cylinder directed opposite the pulsar's velocity vector; 
this region is labeled B2 in Figure^ Material in region 
B 1 generally m oves supersonically , while that in region 
B2 is subsonic (Bucciantini 2002); there thus may be 
significant shear at the interface between these two re- 
gions. Region B is thought to have been observed in ra- 
dio and in X-rays around se veral high-velocity pulsars 
(Chatteriee & Cordes 2002, see also Table 0}, but pre- 
vious observers have generally not made a distinction 
between material in region Bl and that in region B2. 

The shocked pulsar wind materi al is bound ed by a 
contact discontinuity (CD). IChen et all lll996 ) present 
an approximate analytic solution for the shape of the 
CD for the two-layer case appropriate for pulsar bow 
shocks. IWilkinl ( 1996) has derived an exact analytic so- 
lution for the one-layer case, which Bucciantini (2002) 
shows is still a reasonable match to the CDs seen in 
pulsar bow shocks. 

C. Shocked ISM: Beyond the CD, the much denser shocked 

ISM material is advected away from the pulsar, form- 

18 Strictly speaking, a simple rescaling does not allow us to exactly recover 
the situation for Voo = c. However, we expect only small differences between 
a full treatment and the approach adopted here. 



ing a cometary tail bounding the tail containing shocked 
pulsar wind material. This region is in turn bounded by 
a bow shock (BS), at which collisional excitation and 
charg e exchange takes place, generating Ha emission 
fe.gjjones et al 120021 Ic.aensler et al l2002hT). 

These regions and interfaces are all indicated in the simula- 
tion shown in Figure|9] Because Figure ^represents a purely 
hydrodynamic simulation, and does not incorporate the ef- 
fects of a relativistic wind or of magnetic fields, we cannot 
expect exact correspondences between this simulation and our 
data. Nevertheless, by comparison of our high signal-to-noise 
X-ray image with this simulation, we can try to identify in our 
data all the expected components of such a system. 

In particular, we expect the synchrotron emission from a 
bow shock to be sharply bounded on its outer edge by the 
CD. Figure|3]clearly demonstrates that the non-thermal emis- 
sion from G359.23-0.82 shows such a sharp outer bound- 
ary in both the X-ray and radio bands, respectively. Fur- 
thermore, in both radio and in X-rays, this edge shows the 
arc-like morphology expected for the CD seen in Figure [9] 
We therefore identify the eastern edge of the "head" as this 
system's CD, marking the boundary separating the shoc ked 
pulsar wind and the shocked ISM. As estimated in A3. II this 
edge l ies 1 ."0 ± 0"2 from the peak of emission; Cami lo et al.l 
(2002) estimated a similar value from the radio emission 
from this region. We therefore calculate for the Mouse a 
projected radius for the CD in the direction of motion of 
r CD = (0.024 ± 0.005)<i 5 pc. 

4.4. Forward Termination Shock 

We now estimate rf s , the radius of the termination 
shock for ward of the pulsar. Theo r etical expec tations 
are that dvan Buren & McCravl 119881: iBucciantinil l2002t 
Ivan der Swaluw et al.l2003l) : 

r F 

-^-«0.75; (1) 

r cd 

a comparable ratio is seen in Figure [9] We can thus in- 
fer r!p S ps O.OI8C/5 pc, corresponding to an angular separa- 
tion from the pulsar 6j S ss 0"75. Unfortunately, the bright 
X-rays from compact emission in the "head" region (see Fig. 
[8jl prevent us from identifying any features in the image which 
might correspond to this interface. 

Nevertheless, without directly identifying the forward TS, 
we can use our estimate of its radius from Equation Q, com- 
bined with the expectation of pressure balance, to estimate the 
ram pressure produced by the pulsar's motion. Assuming an 
isotropic wind, and in the case where the pulsar's motion is 
wholly in the plane of the sky, we can then write: 

E , 

— rr^pV' W 

47rr£ s c 

where V is the pulsar's space velocity in the reference frame 
of surrounding gas. If the pulsar's motion is inclined to 
the plane of the sky, the situation becomes more compli- 
cated; while the projected separation between the pulsar and 
the apex of the TS will be s maller than the true separa- 
tion dChatteriee & Cordesl2002l) . the projected separation be- 
tween the pulsar and the projected outer edge of the three- 
dimensional surface correspondin g to the TS will be larger 
than the true separation (Gaensler et al. 2002b). Thus in gen- 
eral r F s will always be an upper limit on the true separation 
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between the pulsar and the forward TS, and the ram pressure 
inferred will be a lower limit. We defer detailed simulations 
of this effect for this source, and here assume that all motion 
is in the plane of the sky. 

Since we have estimates of both rj S and E for this sys- 
tem, we can use Equation Q to simply derive that the 
ram pressure produced by the Mouse is pV 2 « 2.1d^ 2 x 
10~ 9 ergs cm" 3 . For cosmic abundances, we can thus write 

V ^ 305n~ l/2 d^ l kms- 1 . 

Rather than assume an ambient density to estimate the ve- 
locity, we can better constrain the properties of the system 
by directly calculating the Mach number, A4. If the speed of 
sound in the ambient medium is c s = V / M. , we can then write: 

pV 2 = M 2 pc 2 = 1ISM M 2 P, (3) 

where -fisM = 5/3 is the adiabatic coefficient of the ISM and 
Pism is the ambient pressure. We adopt a representative ISM 
pressure Pism/^ = 2400Po K cm" 3 , where Pq is a dimen- 
sionless scali ng parameter; typical value s are in the range 
0.5 < Pq < 5 (Feme re1l2()0lllHeilesll2001l) . We thus find that 

the pulsar Mach number is M. = 62c/J 1 P '^ 2 , independent of 
the value of p. This justifies the assumption of a hi gh Mach 
number adopte d in the simulation des cribed in ^4.3l above. 

We note that Yusef-Zadeh & Ballv ( 1987) derived a much 
lower Mach number, M. w 5, by assuming that the edges of 
the broad radio tail seen in Figure [3b) trace the Mach cone 
produced by the pulsar's supersonic motion. However, the 
Mach cone should manifest itself only in the outer bow shock 
jBucciantinil2002l) . which is seen in Ha around some pulsars 
but which is not detected here. 

For the three main phases of the ISM, cold, warm and hot, 
the sound speeds are approximately 1, 10 and 100 km s" 1 , 
respectively. If the pulsar is traveling through cold gas, the 

implied velocity is V = A4c s w 60d^P 1 ' 2 km s" 1 which, 
for 0.5 < Po < 5, i s slower than all but a few percent of the 
overall population ( Arzoumanian et al. 2002). Similarly, if 
the Mouse is embedded in hot gas, the implied velocity is 

V « 6000d^ l P~ 1/2 km s"\ which is well in excess of any 
observed pulsar velocity. We are left to conclude that the 
pulsar is most likely propagating through the warm phase 
of the ISM (with a typical density «o ~ 0.3 cm" 3 ), implying 
a space velocity V w 600c/g 1 P Q km s" 1 . Such a velocity 
is comparable to those for other pulsars with observe d bow 
shocks (see e.g., Table 3 of Chatteriee & Cordes 2002), and 
falls near the center of the expected pu lsar velocity distribu- 
tion at birth l lArzoumanian et ani200 2). The implied proper 

motion is p « 25d^ 2 P 1 ' 2 milliarcsec yr" 1 at a position angle 
(north through east) of w 90°. This is probably too small to 
easily detect with Chandra, but can be tested by multi-epoch 
measurements with the VLA. 

4.5. Backward Termination Shock 

The simulation in Figure |9] shows that the BS and CD are 
both open at the rear of the system, but in contrast, the TS is 
a closed structure. Specifically, the TS is elongated along the 
direction of motion, because while at the apex the pulsar wind 
is tightly confined by the ram pressure of the pulsar's motion, 
in the opposite direction confinement results from pressure 
in the bow shock tail, which can be significantly less. This 
characteristic elongation of the TS is also seen in simulations 
of bow shocks around other systems, such as around runaway 



O stars (Ivan Burenll9 93) and in the Sun's interaction with the 
local ISM <Zankll999ir 

Directly behind the pulsar, the backward termination shock 
is at a distance from the pulsar rf s , where ^> rj S . Thus 
although we have argued above that we cannot observe the 
forward TS, the backward TS might be observable in our 
data. Considering Figure [3ja), we indeed see an elongated 
structure resembling the TS in the simulation, namely the 
"tongue", beyond which the brightness suddenly drops by a 
factor of 2-3 (see Figs. [3 a] and|4). Because the morphol- 
ogy of the "tongue" is very similar to that seen in Figure [9] 
for the cross-section of the TS, it therefore seems reasonable 
that the perimeter of the "tongue" marks the TS in all direc- 
tions around the pulsar. However, there are several important 
differences between the appearance of the "tongue" and that 
expected for the TS in a bow-shock PWN, which we now dis- 
cuss. 

First, the appearance of the "tongue" differs substantially 
from the TS structures seen around the Crab pulsar and other 
young and energetic systems. Most notably, the unshocked 
wind in the Crab Nebula corresponds to a region of mini- 
mal X-ray emission, reflecting the fact that the outflow in this 
area lacks both the energy and the distribution in pitch angle 
to radiate effectively. In contrast, here the region interior to 
the "tongue" is one of the brighter regions of the PWN. Fur- 
thermore, the TS in many Crab-like PWNe appears to tak e 
the form of an inclined torus (see e.g., Ng & Romani 2004), 
suggesting that synchrotron emitting particles are produced 
only in the equatorial plane defined by the pulsar's spin axis. 
(Whether this results from an outflow focused into the equa- 
torial plane, or corresponds to an isotropic outflow for which 
particle acceleration is efficient only in the equator, is a mat- 
ter of debate.) In contrast, the "tongue" does not resemble a 
torus at any orientation, even one that might be elongated by 
the pressure gradient between the forward and backward TS. 

To account for both of these discrepancies, we propose that 
the wind pressure from PSR J1747-2958 is close to isotropic, 
and that particles are accelerated at the TS in all directions 
around the pulsar. In this case, the TS should take the form of 
an ellipsoidal sheath, with the pulsar offset towards one end, 
as shown in Figure |9] No toroidal structure should be seen 
and, because the TS is present in all directions, the underlu- 
minous unshocked wind is completely hidden from view. 

While isotropic emissivity for the TS is not what is 
seen for PWNe such as those around the Crab and 
Vela pulsars or around PSRs B1509-58 and J1811-1925 
JHester et al.ll995t(H elfand et al.l2001tlCiaensler et alJ2 002a: 
Rober tset al.l 120 03). there are various other PWNe imaged 
with Chandra which show a more uniform distribution of 
outflo w and/or i llumination: the PWNe around PSR J 1 1 24— 
5916 (iHughes et al J l200l and PWNe in SNRs G21. 5-0.9 
JSlane et alJl2000t) and 3C 396 dOlbert et al.ll2003l) . all show 
amorphous X-ray morphologies which lack the clear "torus 
plus jets" structure seen for the Crab. The collective prop- 
erties of optical pulsar bow shoc ks also argue for iso tropic 
outflows in those sources (Chatteriee & Cordes 2002). The 
difference between pulsars which show prominent axisym- 
metric termination shocks and others which are surrounded 
by more isotropic structures may be a result of variations in 
the composition of the pulsar wind: regions of efficient parti- 
cle accelera tion may be only tho se in which there are ions in 
the outflow ( Hoshino et ail ll992l) . while the strength of colli- 
mated jets along the spin axis may be a sensitive function of 
the wind magnetization (Komissarov & Lvubarskv 2003). 



X-RAY IMAGING OF THE PULSAR BOW SHOCK G359.23-0.82 



11 



The angular separation between the pulsar and the rear 
edg e of the "tongue" i s &, n n mie . ~ 10", so that 9 t , 
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13. iBuccia ntini (200J) and Ivan der Swaluw et all (120031) ar- 
gue that the wind pressure behind the pulsar is balanced by 
the thermal pressure of the ambient ISM, and correspond- 
ingly provide simple expressions for the ratio of backward 
and forward termination shocks, in their cases giving ratios 
r rs/ r TS A4. However, their formulations are only valid for 
the relatively low Mach numbers used in those simulations. 
We have carried out simulations for a series of higher Mach 
numbers A4 ^> 1, which show that in this regime the pres- 
sure downstream of the pulsar is much higher than that of 
the ISM, and that the ratio of termination shock radii tends 
to an asymptotic limit, r^ s /r^ s ss 5, as can be seen for the 
high Mach number simulation shown in Figure |9] We con- 
clude that the spatial extent of the "tongue" westwards of the 
pulsar, r tongue , is much larger than the expected value of rf s 
— i.e., this region is much more elongated than expected if 
the perimeter of the "tongue" demarcates the TS as proposed 
above. 

Also problematic is that if the "tongue" is a bright hol- 
low sheath, representing the point at which particles are ac- 
celerated at the TS, then we expect it to show significant 
limb-brightening. In contrast, the X-ray emission from the 
"tongue" in FigureQIa) is approximately uniform in its bright- 
ness in the north-south direction, showing only a fading in flux 
as one moves westwards, away from the pulsar (see Fig.@}. 

In the following discussion, we show how both these dis- 
crepancies in the "tongue", namely its excessive elongation 
and lack of limb-brightening, can both be accounted for 
through two effects now observed in many other pulsar winds: 
the finite thickness of the TS, and Doppler beaming of the 
post-shock flow. 

It has been argued that gyrating ions in the TS of a 
pulsar wind can generate magnetosonic waves, which in 
turn can acceler ate electrons an d positrons up to ultrarela- 
tivistic energies (Hoshino et afl [19921) . The resulting shock 
shows significant structure, compression of electrons and 
positrons at the ion turning points resulting in narrow re- 
gions of enhanced synchrotron emission, which may pos- 
sibly account for the "wisps" seen around the Crab Neb- 
ula and around PSR B 1509-58 JGallant & Aronsl IT99H 
Gaens leret alJ200"2al) . In this model, the effective width, Ar, 
of the emitting region at the TS is approximately half of an 
ion gyroradius. We can therefore write: 



Ar = 



1 71 m,c 

2 ZeB 2 



(4) 



where 71 is the upstream Lorentz factor of the flow, ;«,■ is 
the ion mass, Ze is the ion charge, and Bo is the magnetic 
field strength do wnstrea m of the electron TS. In the model of 
Gallant & Arons ( 1994), the upstream Lorentz factor is: 



7i =V- 
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(5) 



where <& ope „ = (E/c) ' is the open field potential of the pulsar, 
and ij is the fraction of this potential picked up by the ions. 
Combining Equations ® and Q, we find that: 



Ar = 



rjE 



?l/2 



2cV 2 B 2 



(6) 



If (Ti is the ratio of energy in electromagnetic fields to that 
in particles in the flow immediately upstream of the electron 



TS dRees & Gunnll974l). th en conservation of energy implies 
JKennel & Coronitl1984bh : 
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where rj$ is the radius of the termination shock in a given 
direction, B\ is the magnetic field just upstream of the electron 
TS, and where in making the final approximation w e have 
assum ed di < 1 (for a review of estimates of o\ . see [Aronsl 
12001 . If (Ti < 1, then B 2 = 3Bi JKennel & Coronitil fl984a). 
Combining Equations and Q, we then find that: 



Ar 

r TS 



6(Tj 



1/2' 



(8) 



Adop ting rj w 1/3 jArons & Tavani|ll993t iGallant & Aronsl 
1994) and <j\ » 0.003 JAronsll2002l) . we then find that the 
thickness of the emitting sheath around the TS should be 
Ar/r rs w 1. 

A key point is that this prediction cannot be maintained at 
the forward TS, where tight confinement of the pulsar wind 
will prevent ions from gyrating; electrons in this region are 
presumably only accelerated at the first ion turning point. 
Thus ahead of the pulsar, we only expect to see emission 
out to a separation r£ s . However in other directions this the- 
ory should be valid, so that the outer edge of the observed 
sheath of emission should be at a separation from the pul- 
sar rjs + Ar « 2rjs- Behind the pulsar, we thus expect a 
maximum angular extent to the "tongue" 9 tongue w 29j S w 
106^ rs 8". This is satisfyingly close to the observed value, 
Otongue = 10", particularly given the approximate nature of the 
above calculation. 

The finite thickness of the emitting region also mitigates the 
effects of limb brightening. For the purposes of the present 
discussion, we crudely approximate the ellipsoidal morphol- 
ogy of the TS by a cylinder, with the thickness of the curved 
walls being equal to the radius of the interior surface. In 
the optically thin case, the factor by which the limbs are 
brightened due to differing path lengths is then: 



fa 



^/(r+Ar^-r 2 
A~r 



V3. 



(9) 



Meanwhile, there is good evidence from PWN morpholo- 
gies that the post-shock flow is sufficientl y relativistic as 
to produce signific ant Doppler boosting (ILu et alJ 120021: 
& Romani 2004). For a downstream flow velocity v 2 = [3c, 
an inclination of the flow direction to the line of sight 4>, and 
an emitting spectrum of p hoton index T, the factor by which 
the flux is boosted is (e.g., Pellin g*et alJl987|) : 



f 



db 



\—(3cos(, 



r+i 



(10) 



In the case a\ <C 1, we have /3 = 1/3 iKennel & Coronitil 
1984a), and we have shown in Tabled that T = 2.1 for the 
"tongue". The elongated geometry of the "tongue" implies 
that along the symmetry axis of the system, the front surface 
of the ellipsoidal emitting region has <f> = 0°, while the back 
surface has <f> = 180°; the northern and southern edges of the 
"tongue" have <fi = 90°. Equation fllOi then implies boosting 
factors of 2.9, 0.34 and 0.83 for the approaching, receding 
and transverse components of the post-shock flow, respec- 
tively. Combining these effects with the difference in path 
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length inferred in Equation (|9j, we find that the observed ra- 
tio in brightness from center to limb in the "tongue" should 
be -1.2. 

In summary, we conclude that the effects known to be oc- 
curring in other PWNe, i.e., finite shock width and Doppler 
boosting of the flow, should produce a sheath of emission in a 
bow shock which has approximately uniform surface bright- 
ness across its extent, and which is twice as elongated as the 
TS seen in simulations, just as is observed for the Mouse. 

4.5.1. Implications for Other Bow Shocks 

Our interpretation that the bright "tongue" of G359.23-0.82 
marks the surface of the TS can also be applied to other bow- 
shock systems — provided that the Mach number is high, we 
expect to similarly see an elongated TS region like that seen 
here. 

In Tabled we list estimates of the observables rcD and l la u 
for all confirmed X-ray bow shocks, where l m a is the maxi- 
mum extent of the X-ray tail in the direction opposite that of 
the pulsar's motion. What is striking from this Table is that 
the ratio Itaii/rcD is significantly larger for the Mouse than 
for other systems. For most of these other pulsars, the ratio 
hail I r cd is comparable to the ratio r tongue /rcD ~ 10 seen here 
(where r tongue is the spatial extent of the "tongue" westwards 
of the pulsar). 

We therefore propose that the elongated X-ray trails seen 
extending opposite the pulsars' direction of motion are not 
necessarily shocked particles well downstream of the flow, as 
is usually interpreted, but rather in some cases m ight repre- 
sent the surface of the TS. (We note that Gvaramadzel l2004l 
has recently come independently to the same conclusion in 
the case of PSR B1757-24.) It is only in a source such as 
the Mouse, with its much higher X-ray flux, that the fainter 
region of downstream flow (the "tail" in Fi gure |3ta] or region 
B in Figure |9j to be discussed further in q4.7l below) can be 
observed. 

We can immediately make specific predictions as to what 
properties other systems should have, depending on whether 
we are seeing the surface of the TS or the flow of postshock 
material downstream of the TS. 

If the emission which we see represents the surface of the 
TS, then we expect that the X-ray spectrum should show lit- 
tle variation across the extent of the source, and since radio 
and X-ray emit ting particles are both apparently accelerate d 
in this region JBietenholz 'et al.ll2001t Idaensler et alJ l2002a). 
radio and X-ray emission should both show a similarly nar- 
row and elongated morphology. In particular, there should 
be an abrupt turn-off or sudden decrease in radio and X-ray 
emission in the direction opposite the pulsar's motion, at an 
extent corresponding to l ta u / rcD ~ 10, marking the position of 
the backward TS. We lack sufficient counts in most of these 
systems to constrain any spectral variation as a function of 
position, but certainly PSR B 1957+20, PSR B 1757-24 and 
CXOU J06 1705. 3+222 127 in TableHappear to meet most of 
the other criteria. We therefore propose that in these sources, 
the cometary emission previously identified corresponds sim- 
ply to the surface of the TS, as is seen in the "tongue" of 
the Mouse. (We note that in the case of PSR B 1757-24, 
iKaspi et al.l2001l use the length of the X-ray trail to calculate 
the flow velocity in shocked wind material. However, in our 
interpretation, all points in this X-ray trail are regions of fresh 
particle acceleration, and such data cannot be used to calcu- 
late flow speeds.) Deeper observations of these systems can 
potentially reveal X-ray emission at greater extents, which in 



our interpretation would correspond to the post-shock flow. 

On the other hand, if the synchrotron emission from a pul- 
sar bow shock corresponds to shocked material downstream, 
then because of increasingly severe but energy dependent syn- 
chrotron losses further downstream, we expect that the X- 
ray spectrum should gradually soften with increasing distance 
from the pulsar, that the radio extent be significantly longer 
than the X-ray extent, and that the X-ray and radio "tails" 
should gradually fade into the background with increasing 
distance from the pulsar. We also expect that the X-ray mor- 
phology should be elongated such that l, a a / rcD 3> 10, and that 
the radio emission should fill a broad region filling the space 
between the TS and CD, corresponding to region B 1 of Fig- 
ure|9] 19 The sources in Table|4]which appear to meet all these 
criteria are the Mouse, PSR B1951+32 and PSR B1853+01. 
For these three sources, we therefore argue that the furthest 
extent of X-ray and radio emission behind the pulsar represent 
postshock material flowing downstream. For PSR B 195 1+32, 
no detailed spectral information is yet available, but a recent 
Chandra image (Moon et al. 2004) indeed suggests the pres- 
ence of two components to the X-ray emission behind the 
pulsar, with similar morphologies to the "tongue" and "tail" 
seen in X- r ays fo r the Mouse. The X-ray observations of 
iPetre et all (120021) lack the spatial resolution and sensitivity 
to constrain the possible presence of a "tongue" of brighter 
emission closer to PSR B1853+01, possibly representing the 
surface of the TS as seen here for the Mouse. Deep on-axis 
Chandra observations are required to investigate this possibil- 
ity. 

4.6. The Halo 

The X-ray image and one-dimensional profile in Fig- 
ures |3 a) and 0] respectively, both clearly show that there is 
a "halo" of emission extending several arcsec eastwards of 
t he s harp drop-off in brightness which we have identified in 
34.3l with the CD. Comparison with Figure|9]suggests that the 
X-ray emission from the halo might correspond to part of re- 
gion C of the simulation, in which ambient gas is heated and 
compressed by the bow shock. We here consider whether the 
X-ray emission seen from the "halo" could be from shock- 
heated material in this region. 

The rate at which kinetic energy associated with the pulsar's 
motion is converted into thermal energy behind the bow shock 
is: 

L h0W = 7rr 2 CD pV^^E, (11) 

where we have used Equations Q and in deriving the fi- 
nal expression. For V = 305n^ 2 d^ 1 km s" 1 as derived in 34.41 
we thus find L/, OM ss lAn^^d^ 1 x 10 33 ergs s" 1 , implying an 

-1 12 

isotropic unabsorbed flux at Earth of fb OK ~ 3.7n d 5 x 
10~ 13 ergs cm -2 s _1 . This value is broadly consistent with the 
flux inferred for region 6 from thermal models with Nh fixed 
at 2.7 x 10 22 cm" 2 , as listed in Tabled Thus from an energet- 
ics argument alone, the expected luminosity of shock-heated 
gas in this region is consistent with what is observed in the 
X-ray band. 

However, there are several reasons why we think it unlikely 
that the emission seen here comes from heated ISM material. 
First, the best spectral fits in TableOJimply an absorbing col- 
umn which is inconsistent with that seen for all of regions of 

19 X-ray emission will generally not show this broad a ppear ance due to 
synchrotron losses in this region; see detailed discussion in 34.71 
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1-5. When one fixes Nh = 2.7 x 10 22 cm -2 to match regions 
1-5, the resulting spectral fits for region 6 are quite poor. Sec- 
ond, the expected gas temperature behind a shock of this ve- 
locity is kT w Q.2n^ l d^ 2 keV. For no = 0.3 cm -3 and d$ = 1 
as adopted here, we infer kT « 0.6 keV, well below the value 
kT > 4 keV seen in thermal fits to region 6 as listed in Tabled 
Finally, even if we have overestimated the ambient density 
and pulsar velocity so that the expected shock temperature is 
indeed kT 4 keV, the expected X-ray luminosity in shocked 
gas near the head is L x = An^V w 4ngrf| x 10 27 ergs s" 1 , where 
A w 10~ 23 ergs cm 3 s" 1 is the cooling function at kT w 4 keV, 
and V ~ 71-r 3 ^ = 4.1d| x 10 cm 3 is the approximate volume 
of shocked gas in the nose of the bow shock. The implied X- 
ray flux at Earth is many orders of magnitude below what is 
observed in this region. 

We rather propose that the emission seen in the "halo" 
corresponds to a combination of dust-scattering and contri- 
butions in this region due to the wings of the PSF For an 
unresolved source with a spectrum as given for region 1 in 
Table [2 one expects the dust-scattered flux falling into re- 
gion 6 to comprise ~ 0.7% of the incident flux, if dust is 
distributed homogeneously between the source and observer 
( R. K. S mith, 2003, private communication, using the method 
of Mathis & Lee 1991). We have used the PSF generated by 
CHaRT (see ^3~4l above) to estimate that — 0.9% of the inci- 
dent flux from a compact source with the spectrum of region 1 
should fall into region 6 resulting from the wings of the PSF. 

Since region 1 contains 7200 counts, we thus expect ~ 120 
counts in region 6 due to these two effects. Although this is a 
factor of ~ 3 less than the cou nt-rate actually observed in re- 
gion 6, the discussion in §4.1l suggests that it is highly likely 
that the scattering medium is located in the molecular ring 
rather than is distributed evenly along the line of sight, and is 
thus concentrated much closer to the source than to the ob- 
server. Incorporating this effect would significantly increase 
the estimate of the contribution of dust scattered emission in 
region 6. Furthermore, our estimate does not take into account 
that region 1 is heavily piled-up, or that there is other bright 
extended emission in the immediate vicinity of region 1, the 
photons from which will also be scattered by dust and by the 
Chandra mirrors. We thus conclude that the 120 counts es- 
timated above is likely a significant underestimate, and that 
the flux in the "halo" region is indeed broadly consistent with 
being light scattered by both dust and by the Chandra mirrors. 

We note that because the cross-section for dust scattering 
decreases with increasing energy, one expects the spectrum of 
dust-scattered emission to be softer than that of the source. 
Indeed for the power law fits to region 6 in Table [2 the ab- 
sorbing column is less or the photon index is steeper than for 
region 1, either result indicating a softer spectrum as expected. 

4.7. The Tail 

The "tail" is perhaps the Mouse's most puzzling aspect. As 
seen in Figures^and|3lb), at radio wavelengths the tail flares 
considerably about 10" west of the pulsar, then further down- 
stream contracts, broadens, and then contracts again into a 
narrow cylindrical structure, 12' in length, which eventually 
fades into the background. Here we attempt to explain the be- 
havior close to the "head", and the long thin radio tail. We do 
not attempt to account for the complicated radio morphology 
between these two regimes, which perhaps reflects inhomo- 
geneities in the surrounding medium (e.g.. Uones et alJ 2002). 

We expect that the X-ray and radio emission from the "tail" 
should correspond to region B of Figure [9] where shocked 



pulsar wind material emits synchrotron emission and flows 
out behind the pulsar. With increasing distance from the pul- 
sar, we expect that the X-ray emission should fade and that its 
photon index should steepen, both as a result of synchrotron 
losses. These effects are both clearly observed in Chandra 
data. However, the properties of the "tail" differ from those 
expected, in that the X-ray "tail" is reasonably narrow and 
collimated, in contrast to the increasingly broad region of 
shocked material seen in region B of Figure |5] 

In the radio image shown in Figure |3jb), the outer regions 
of the tail do appear to flare and broaden, showing behavior 
broadly similar to t hat seen in the simulation. However, as 
pointed out in 33. 21 and shown in Figure |5] the radio tail con- 
sists of two components: the broad region just discussed, plus 
a brighter narrower component closer to the symmetry axis, 
and showing very good morphological correspondence with 
the X-ray image. 

Comparison with Figure [9] suggests that we indeed expect 
two components to the "tail": region Bl, in which shocked 
material flows from the forward TS, and region B2, where 
material shocked at the backward TS flows away from the 
pulsar in a narrower region. We propose that the outer and 
inner components of radio emission from the "tail" seen in 
Figures |3jb) and [5] correspond to regions Bl and B2 of the 
simulation, respectively, and that only region B2 is seen in 
X-rays. 

For any reasonable value of the nebular magnetic field, the 
synchrotron lifetime of X-ray-emitting electrons will be sig- 
nificantly less than the age of the system. The X-ray extents 
seen in the post-shock flow thus act as a useful diagnostic of 
the flow speed and magnetic field strength in each region. 

Although our hydrodynamic simulation does not explicitly 
include relativistic effects or magnetic fields, both of which 
are likely to play an important role in moderating the proper- 
ties of the post-shock flow, we can use the simulation to make 
a rough estimate of the relative flow velocities in each of re- 
gions Bl and B2. We find in the simulation that material in 
region Bl has an asymptotic bulk velocity ~ 4 times larger 
than in region B2. 

This difference in velocities can be understood by the fact 
that most particles feeding region Bl do not cross the termina- 
tion shock head-on, and are accelerated by pressure gradients 
in the head of the bow shock. In contrast, particles feeding 
region B2 meet the shock head on, propagate through an es- 
sentially uniform pressure distribution, and so are not accel- 
erated. Thus if the magnetic field strengths were the same in 
regions B 1 and B2, we would expect X-ray emission to extend 
far from the pulsar in a broad region corresponding to B 1 , and 
would see a much shorter scale for the narrow component of 
X-rays corresponding to B2. However, we clearly observe the 
opposite situation in this system: despite corresponding to a 
slower bulk flow speed, the spatial extent of X-rays in region 
B2 is significantly larger than for Bl. 

A simple interpretation is that the magnetic field strength 
in the post-shock flow is markedly lower in region B2 than 
in region Bl, which then allows X-ray-emitting electrons to 
propagate much further from the shock before radiating away 
their energy. We quantify this as follows. We first estimate 
the length scales of the X-ray emitting regions. We estimate 
that in region B 1 , X-ray emitting material extends about half 
way around the northern and southern sides of the "tongue", 
but is not detectable further downstream. The corresponding 
angular distance traveled by X-ray emitting material is w 8", 
or a distance of x F w Q.2ds pc. For region B2, we estimate 
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from Figures and |4] an angular extent ~ 40" in X-rays, 
corresponding to a spatial extent x B w l.Ods pc. 

If the length scale of X-rays in each region is determined 
by the time traveled before synchrotron losses become signif- 
icant, then 

^!!flf (12) 

x B V B \B B J ' K } 

where V is the flow velocity, B n is the nebular magnetic field, 
the superscripts F and B correspond to situations for par- 
ticles flowing from the forward and backward termination 
shocks respectively. From the above discussion, we have that 
x F Jx B w 1 15 and V F /V B « 4; Equation then implies that 

Since our simulation does not explicitly include magnetic 
fields, we do not attempt to calculate actual values of the field 
strength in each region. However, we note that if the degree 
of magnetization of the post-shock flow is similar in the for- 
ward and backward directions, then an expression similar to 
Equation Q implies that 

R F B 

^L= r -1± (13) 

B B r F ■ ^-^ 
D n r TS 

Using rj S m 0.018^5 pc and r B s w 0.5r tongue = 0.12d 5 pc, we 
then find B F /B B as 7, in good agreement with the estimate 
made above from the X-ray morphology. While this simple 
treatment provides support for the argument that the X-ray 
morphology of the tail is governed by differing magnetic field 
strengths in the different parts of the tail, a full relativistic 
magnetohydrodynamic treatment of this system is needed to 
properly address the amplification and evolution of magnetic 
fields in each region. Such an analysis is beyond the scope of 
this paper, but will be the subject of a future study. 

If synchrotron losses indeed limit the extent of the X-ray 
"tail", we expect the tail to be increasingly shorter in extent at 
progressively higher energies. We have formed images of the 
Mouse in various energy bands (0.5-0.8, 0.8-2.0, 2.0-5.0 and 
5.0-8.0 keV), and qualitatively can confirm that indeed the tail 
appears to shrink as one considers successively higher energy 
photons. However, we lack sufficient counts to quantify this 
effect, as the energies where we expect the tail to be longest 
are those where interstellar absorption is most severe. 

Synchrotron lifetimes are proportional to e" 1 / 2 , where e is 
the energy at which photons are emitted. We thus expect the 
lifetime of radio-emitting particles in regions B 1 and B2 to be 
~ 10 4 times longer than their X-ray counterparts. We thus ex- 
pect negligible synchrotron losses in the radio band. Through- 
out region Bl, we thus expect to see considerable emission 
at radio wavelengths. Indeed this expectation is fulfilled — 
while there is little X-ray emission in much of region B 1 , there 
is clearly significant radio emission in the increasingly broad 
component of the "tail" seen far from the symmetry axis in 
Figure |3 Further downstream, adiabatic losses may eventu- 
ally become significant in the radio-emitting particles as the 
flow diverges, while changes in the density of the ambient 
medium may have a dramatic effect on the morphology of 
this region. 

To consider the properties of radio emission in region B2, 
we assume a post-shock velocity V B = c/6 at the backward ter- 
mination shock, in the reference frame of the pulsar (L vutikovl 
2004). Downstream of the backward termination shock, the 
flow velocity evolves as V B oc 1 /A, where A is the cross- 
sectional area of region B2. The morphology of radio emis- 
sion in Figures ^and|3lb) suggests that the width of the radio 



"tail" is a factor of ~ 3 broader than the width of the back- 
ward termination shock as defined by the morphology of the 
X-ray "tongue". We therefore adopt a constant flow veloc- 
ity V B f» c/6/3 2 ~ 5600 km s" 1 in the far downstream com- 
ponents of region B2, so that in an external reference frame 

the flow moves westwards at a speed (5600- 600^5 'Pq 1//2 ) ~ 
5000 km s" 1 . 

We have not included magnetic fields in our simulation. 
However, we expect that magnetic field stresses should main- 
tain collimation of the flow far downstrea m in region B2, 
as is seen for the Earth' s magnetotail (e.g., Bavassan o et alJ 
119741 ISlavin et alJll9 83) and for simulation s of old n eutron 
stars moving through the ISM ( Toropina et al. 200 0). In 
this case adiabatic, as well as synchrotron, losses in the ra- 
dio band will be negligible, and emission from the tail can 
then have a consider a ble ex tent. Indeed the radio data of 
Yusef-Zadeh & Ballv (1987) show that the tail extends for 
ps 12' west of the pulsar, corresponding to an extent x 8 sa 

lids pc. At a pr ojected velocity V w 600t/J 1 f 1 ^ 2 km s" 1 as 
inferred in SI4.4I and assuming that the pulsar's true age is 
t = 25.5 kyr, the total projected distance traveled by the pul- 
sar in its lifetime is L sa I6d$ l P 1//2 pc. Since L w x 8 , it seems 
reasonable to suppose that the pulsar was born near the tail's 
western tip, with the tail representing a wake of particles left 
behind by the pulsar as it travels. Such a possibil ity has been 
suggested previously by other authors JPredehl & Kulkarnll 
H995llCamilo etal 120021) . 

However, this suggestion is problematic for several reasons. 
First, the radio data in Figure^indicate that there is no sud- 
den termination of the tail at its westernmost extent. Rather, 
beyond this point the tail simply becomes too faint to see, 
particularly in the presence of emission from SNR G359.1- 
0.5. Second, the calculation made in the preceding para- 
graph neglects the backward flow velocity of emitting parti- 
cles, which should cause the tail to extend further westwards 
of the pulsar's birthplace as a function of time. If we adopt 
a constant flow velocity V B « 5000 km s" 1 in region B2 as 
calculated above, the full length of the bow shock structure 
could be many times longer than its observable component, in 
which case the fact that L sa x 8 would simply be coincidence. 
Third, although no associated SNR has yet been identified for 
PSR J1747-2958 (perhaps because the SNR ex panded into a 
low density environment or pre-existing cavity; Kafat oTet alJ 
119801 iBraun et"afl 11989ft . the supernova blast wave should 
have a significant impact on its environment, sweeping up 
surrounding material and evacuating a large cavity around the 
pulsar birthplace. If the pulsar was indeed born at the tip of 
the radio tail of the Mouse, it is hard to understand how the 
tail extends smoothly and uninterrupted through the turbulent 
and highly inhomogeneous environment created by the asso- 
ciated SNR. Fin ally, for V sa 600 km s" 1 and no sa 0.3 cm -3 as 
adopted in 34.41 the pulsar should have broken out of its asso- 
ciated SNR when the latter wa s still in the Sedov-Taylor p hase 
of evolution (see Eq. (4) of Ivan der Swaluw et al]l2"003l) . In 
this case, it is simple to show that the age of the pulsar when 
it escapes from its SNR is: 

tescape = UE^ (V / 1000 km s" 1 r 5/3 «~ 1/3 kyr, (14) 

where £ V| x 10 51 ergs s" 1 is the kinetic energy of the ex- 
plosion Jshull et alJH989t Ivan der Swaluw et alJl200l . For 

V sa 600iiJ 1 f > 5 ^ 6 km s" 1 and no sa 0.3 cm" 3 , we find t escape m 

49E l 5 { 3 d 5 5 /3 P ~ 5/6 kyr. If the true age of PSR J1747-2958 is 
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t = 25.5 kyr, then for £51 ~ 1, ^ « 1 and Pq w 1, the pul- 
sar should still be embedded in its associated SNR, which 
is not the case. A simple resolution for this discrepancy is 
if PSR J1747-2958 is significantly older than its character- 
istic age, as has also been claimed for P SRs B0833^-5 and 
B 1757-24 llLvne et alJfT9%t iGaensler & Fraill EOQO). The 
pulsar birthplace would then be significantly further west of 
the endpoint of the radio tail. 

To summarize, the radio and X-ray morphologies of the tail 
can both be consistently explained by a model in which there 
are two regimes of downstream flow, corresponding to mate- 
rial shocked at the forward and backward termination shocks. 
At the forward TS, a high magnetic field causes severe syn- 
chrotron losses, so that X-ray emitting particles in this region 
lose their energy before traveling only a fraction of a parsec; 
radio-emitting particles form a broad cometary region, match- 
ing the morphology of this region seen in simulations. At the 
backward TS, the flow speed is somewhat lower than at the 
forward TS, but the magnetic field is very weak. The result- 
ing synchrotron lifetime and the minimal divergence of the 
flow allow a collimated and significantly extended X-ray tail 
to form. At radio wavelengths, the flow from the backward 
TS is what produces the long narrow tail for which the Mouse 
originally attracted attention. The visible endpoint of this tail 
probably does not mark the pulsar birthplace; most likely the 
pulsar is much older than its characteristic age, and was born 
well beyond the point where radio emission from the tail be- 
gins. 

5. CONCLUSIONS 

We have here presented a detailed study of the X-ray emis- 
sion from G359.23-0.82 ("the Mouse"), along with an im- 
proved timing solution for the coincident young radio pulsar 
J1747-2958. Our main findings are as follows: 

1. G359.23-0.82 is a bright X-ray source with a 
cometary morphology and a power-law spectrum, with 
PSR J1747-2958 embedded in its head. This conclu- 
sively demonstrates that the Mouse is a bow shock pow- 
ered by the relativistic wind of PSR J1747-2958, all at 
a likely distance of w 5 kpc. The Mouse is a highly 
efficient X-ray emitter, converting 2% of its spin-down 
luminosity into synchrotron emission at these energies. 

2. High-energy emission from the pulsar itself is likely 
identified in the form of a compact X-ray source co- 
incident with the pulsar position. This most probably 
corresponds to pulsed emission from the pulsar magne- 
tosphere. A slightly extended clump of uncertain nature 
lies immediately adjacent to the pulsar, and produces 
> 30% of the total X-ray flux from this system. 

3. Using the high spatial resolution of Chandra data, we 
identify the contact discontinuity separating shocked 
wind material from shocked ISM ahead of the pulsar's 
motion. Pressure balance then allows us to infer that 
the pulsar's Mach number is M. « 60, implying a likely 
space velocity V ~ 600 km s" 1 through warm gas with 
density «o ~ 0.3 cm" 3 . 
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4. We see a bright elongated region extending 10" behind 
the pulsar, which we identify as a sheath of bright X- 
ray and radio emission surrounding the wind termina- 
tion shock of an isotropic flow. A model which incor- 
porates both the finite thickness of an ion-loaded wind 
shock and Doppler beaming of the relativistic flow can 
account for both the size of this region and its lack of 
limb brightening. We argue that the X-ray and radio 
trails seen behind PSR B 1957+20, PSR B 1757-24 and 
CXOU J06 1705. 3+222 127 similarly represent the sur- 
face of the wind termination shock in these systems. 

5. The tail of the Mouse contains two components: a col- 
limated region seen in both X-rays and in radio, and a 
broader component seen in radio only. We propose that 
the properties of the former are consistent with emis- 
sion produced by the post-shock flow beyond the back- 
ward termination shock behind the pulsar, where the 
magnetic field is low and synchrotron losses are mod- 
erate. The latter likely corresponds to material shocked 
in a high-field region ahead of the pulsar, where radia- 
tive losses are severe at X-ray energies. Similar behav- 
ior is possibly also seen in the radio/X-ray tails behind 
PSRs B1951+32 and B1853+01. 

Taking advantage of the high X-ray count-rate seen from 
this source, we have successfully used the Mouse as a testing 
ground for combining our understanding of the shock struc- 
tures seen in both Crab-like and bow-shock PWNe. How- 
ever, we caution that while recent Chandra observations have 
showed some satisfying similarities among PWN properties, 
they also reveal some puzzling differences. It is clear that to 
properly understand the pulsar bow shock phenomenon, we 
will need to obtain deep observations of other such systems, 
which can then be compared and contrasted to the remarkable 
Mouse. 
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TABLE 1 

Spectral fits to regions 1-6 of G359.23-0.82. 



Region 


Description 


Total Counts 
(0.5-8.0 keV) 


Nh 
(10 22 cm 4 ) 


F/kT 
— /keV 


fx 

(10 12 ergs cm 2 s ') a 


2 / 


1 

2 
3 
4 
5 


Head (inner) 
Head (outer) 

Tongue 
Tail (inner) 
Tail (outer) 


7237 ±85 b 
1041 ±32 
1457 ±38 
1100 ±33 
1143 ±34 


2.7 ±0.1 


1.8±0.1 
1.9±0.1 
2.1±0.1 
2.3±0.1 
2.5±0.1 


11.2±6.6 b 

1 2 +03 
1 - -0.2 

1.8 ±0.3 
1.5 ±0.3 

1 -'-0.3 


361/312 = 1.16 


6 
6 
6 
6 


Halo (PL) 
Halo (RS) 
Halo (PL) 
Halo (RS) 


315 ±19 


Lo -0.5 

1.9±0.5 
2.7 (fixed) 
2.7 (fixed) 


1 5 +l, ° 
'■ J -0.4 

>4.4 
2.1 ±0.3 
4 1 +1J 


<0.06 
0.03 ±0.01 
0.4±0.1 
0.33 ±0.05 


8.1/8 = 1.0 
6.7/8 = 0.8 
12.1/9= 1.4 
12.7/9= 1.4 


Note. 


— Uncertainties 


are all at 90% confidence. All models assume interstellar absorption 


using the cross- 



sections of Baiuciriska-Church & McCammon 1 1992), assuming solar abundances. Regions 1-5 have been fit with 
a power law model, of the form f £ <x e~ r where Y is the photon index. Regi on 6 is fit w ith both a power law 
("PL") and a Raymond-Smith spectrum of temperature T ("RS") i Raymond & Smith 1977). The model fitted to 
the spectrum of region 1 also includes the effects of pile-up, according to the prescription of Davis 1 2003). 



a Flux densities are for the energy range 0.5-8.0 keV, and have been corrected for interstellar absorption. 
b The number of counts quoted for region 1 is that detected, i.e., this value suffers from pile-up. The flux density 
quoted for region 1 is that incident on the detector, i.e., this value has been corrected for the effects of pile-up. 



TABLE 2 

Best-fit parameters for the thre e com ponent model for the "head" 
ofG359.23-0.82,as discussed in ii3.4l All uncertainties represent 
90% confidence intervals. 



Model Component 


Parameter 


Best-Fit Value 


1st Gaussian 


RA of Center (J2000) a 
Deck of Center (J2000) 3 
FWHM (arcsec) 
Ellipticity b 
Position Angle 

of Major Axis (N through E) 
Peak Amplitude (counts arcsec -2 ) 


17 h 47 m 15?854± 0*004 
-29°58'0i;'38±0;'04 
1.1 ± 0.1 
0.25 ±0.09 
118° ±12° 

1320 ± 150 


2nd Gaussian 


RA of Center (J2000)" 
Deck of Center (J2000f 
FWHM (arcsec) 
Ellipticity b 
Position Angle 

of Major Axis (N through E) 
Peak Amplitude (counts arcsec -2 ) 


17"47 m 15?765± 0*005 

-29°58'01"29±0?'04 

2.4±0.1 

0.12 ±0.06 

91° ± 13° 

925 ±73 


Offset 


Amplitude (counts arcsec -2 ) 


10± 13 



a The unc ertai nties in position are those within the X-ray reference frame. As dis- 
cussed in j|3.1l absolute positional uncertainties in this frame are fa 0"25 in each co- 
ordinate. 

b In SHERPA, the ellipticity, e, is defined by e = l — b/a, where a and b are the widths 
along the major and minor axes respectively. 
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TABLE 3 

Pulsar parameters for PSR J 1 747-2958 , as 

MEASURED AND INFERRED FROM RADIO TIMING 
OBSERVATIONS. 



Parameter 


Value 


RA (J2000) 


17 h 47 m 15?882(8) 


ueci. (jzuuuj 


—Zy Jo Ul. \)\ / ) 


Rotation frequency, v (s -1 ) 


10.1200258870(2) 


Frequency derivative, v (s -2 ) 


-6.28042(1) x 10 -12 


Frequency second derivative, i) (s -3 ) 


9.1(2) x 10~ 23 


Epoch (MJD [TDB]) 


52613.0 


Data span (MJD) 


52306-52918 


R.M.S. residual (ms) (white/red) . . . 


0.5/4.8 


Dispersion measure, DM (cnT 3 pc) 


101.5 


Derived parameters: 




Characteristic age, r (kyr) 


25.5 


Spin-down luminosity, E (ergs s~') 


2.5 x 10 36 


Magnetic field strength, B (G) .... 


2.5 x 10 12 



NOTE. — Numbers in parentheses represent the uncertain- 
ties in the last digits quoted and are equ al to the formal 1-cr 
errors determined with TEMPO. See i|3.5l for an explanation of 
how these parameters were obtained. 



TABLE 4 

The known sample of X-ray bow shocks around pulsars. 



Source 


Other Name Distance 


Projected Velocity 11 


r C D b 


hail 


kai/rcD L x (0.5-8.0 keV) 


Lx/E 


Reference 




or SNR (kpc) 


(km s" 1 ) 


(pc) 


(pc) 


(ergs s~') 







PSR J 1747-2958 


Mouse 


5ch 




0.024^5 


> 1.1* 


>46 


5rf| X 10 34 


2dj x ur 2 


PSR B 1951+32 


CTB 80 


2 


240 ±40 


0.02 


"^0.4 


~20 


~2x 10 33 


~ 5 x 10~ 4 


PSR B 1853+01 


W44 


2.6 




~ 0.05 


~ 1 


-20 


~ 5 X 10 32 


~ 1 x 10~ 3 


PSR B 1757-24 


Duck 


5 


< 590 


0.036 


0.5 


14 


~ 4 x 10 32 


~1.4x 10"" 4 


PSR B 1957+20 


Black Widow 


1.5 


220 ±40 


< 0.01 


0.12 


> 12 


~ 3 x 10 30 


- 3 x 10- 5 


CXOU J061705.3+222127 


IC443 


1.5 




0.06 


~0.5 


~8 


~5x 10 33 





This paper 
Moon et al. 2004 
Petre et al. 2002 
Kaspi et al. 2001 
Stappers et al. 2003 
Olbert et al. 2001 



Note. — We are aware that there are several other pulsars for which X-ray bow shocks have been claimed in the literature, but here list only those symciiis whose idemhicalion as a how shock seems secure. Of X-ray bow shocks listed in Table 3 of 
Chatter jce & Cordcs 2002 i and elsewhere, we exclude here: PSR B2224+65. for which the original faini X-ray detection was not seen in subsequent observations Wong ci al 2002 I; PSR B 1823-1 3, for which more sensitive observations have shown that no 
bow-shock morphology is present (Jaenslcr et al 2003 I; PSR B 1929+10, for which the Taint X-ray trail seen by ROSAT is not apparent in deeper observations (M. J. Pivovaroff. 2003. private communication): and PSR 10537-6910. for which a bow-shock 
interpretation is problematic Irouler^raE^BBwnVOT^er^waluw etaHZWM. We have also omitted various other sources, including PSR J1016-5857 ICamilo et al|200H . G0.13-0.ll I Wane et ald2BBl . GeV J1809-2327 IBraieetallSoT^ and Geminga 



Caraveo el al 2003 i. all of which have been postulated as possible X-ray bow shocks but whose nature has not yet been confirmed. 

a Projected velocities are only listed in cases where they have been measured directly from proper motion or timing measurements. 

Radius of the contact discontinuity ahead of the pulsar. 
C Maximum X-ray extent in the direction opposite the measured/inferred direction of motion. 
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17 h 47 m 30 s 15 s 47 m s 45 s 30 s 15 s 46 m s 

RA (J2000) 

FIG. 1. — Existing radio and X-ray data on G359.23-0.82 ("the Mouse") and PSR J1747-2958. The image shows a Very Large Array (VLA) image at a 
frequency of 1.5 GHz. The brightness scale is logarithmic, ranging between -1.3 and 87.6 mJy beam -1 as shown by the scale bar to the right of the image. This 
image is made from VLA observations in the CnB and DnC configurations, carried out on 1986 Oct 11/12 and 1987 Feb 20, respectively. The spatial resolution 
is 12" 8 X 8"4, a s shown by the filled ellipse in the lower right corner of the image. The original l-cr error ellipse on the position of PSR J1747-2958, as reported 
by Camilo et al. 1 2002), is shown at the head of the radio nebula. Also overlaid are contours from a ROSAT PSPC observation carried out on 1992 Mar 02 of 
duration 2000 s, convolved with a gaussian of dimensions 20" X 20". Contour levels are drawn at levels of 4%, 8%, 20% and 60% of the peak. The two straight 
lines indicate the southern and western edges of the Chandra ACIS-S3 field of view. 
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17 h 47 m 36 s 30 s 24 s 18 s 12 s 6 s 47 m s 54 s 



RA (J2000) 



FIG. 2. — Chandra image of the ACIS-S3 field of view in the energy range 0.5 to 8.0 keV. This image has not been corrected for spatially-varying exposure 
across the field of view. Sources discussed in the text are labeled. The region lying between the two circles (but excluding G359.23-0.82, the USNO-A2.0 star 
0600-28725346 and the read-out streak from SLX 1744-299) was used to estimate the background level when carrying out spectroscopy of G359.23-0.82. 
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FIG. 3 . — X-ray and radio images of the Mouse, (a) Chandra image of the Mouse in the energy range 0.5-8.0 keV. The brightness scale is logarithmic, ranging 
between and 100 counts, as shown by the scale bar to the right of the image. Various features discussed in the text are indicated. This image has not been 
exposure corrected, (b) VLA image of the Mouse at a frequency of 4.8 GHz. The brightness scale is linear, ranging between -0.2 and 6.0 mJy beam -1 as shown 
by the scale bar to the right of the image. The contours are at levels of 0.9 and 1.8 mJy beam -1 , chosen to highlight the two different components seen in the 
radio "tail". This image is made from VLA observations in BnA and CnB arrays, carried out on 1987 Nov 07 and 1988 Feb 25, respectively, with a resulting 
spatial resolution of 1"6 X 1"7. 
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1000 F 1 1 1 1 1 1 1 1 1 1 r 




Offset (arcsec) 

FIG. 4. — 0.5-8.0 keV brightness profile of G359.23-0.82 in the east-west direction (negative offsets are to the east). The solid line and data points represent the 
X-ray surface brightness of G359.23-0.82 as a function of offset from the peak, averaged over a column of nine 0"492 X Of'492 pixels bisected by the symmetry 
axis of the source. The dashed line indicates the brightness of the background, determined by averaging a column of 121 pixels immediately to the north and 
south of G359.23-0.82. The dot-dashed line divides the logarithmic scale in the upper panel from the linear scale in the lower panel. Various features discussed 
in the text are indicated. 
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Offset (arcsec) 

FIG. 5. — Radio and X-ray brightness profiles of the "tail" of G359.23-0.82 in the north-south direction (negative offsets are to the south). The solid line 
represents the X-ray surface brightness of the tail as a function of offset from the symmetry axis, averaged over a row of 26 Of'492 X ()"492 pixels centered 
on RA (J2000) 17 h 47 m 13?8. The dashed line indicates the radio surface brightness of the tail as a function of offset from the symmetry axis at RA (J2000) 
17 h 47 m 13?4. Both profiles have been Hanning smoothed with a smoothing length of five ACIS pixels, or 2"46. The narrow component of the tail is seen in both 
radio and X-rays, while the broad component appears only in the radio. 
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FIG. 6. — Image as in Fig.|3ja), but showing the extraction regions used for spectroscopy. The five regions indicated are mutually exclusive, each lying between 
successive sets of contours, but consisting only of data lying to the west of the vertical line. Region 6 is not shown, but consists of all points to the east of this 
line which fall within an annulus centered on the peak, of inner radius 4" 8 and outer radius 14" 1. 
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0.5 



Channel Energy (keV) 



FIG. 7. — Chandra spectrum of G359.23-0.82 in the five different regions shown in Fig.|6] The points in the upper panel indicate the data in each region, while 
the solid lines show the corresponding best fit power law models, resulting from fits to the data in which all regions are constrained to have the same value of Nfj . 
The spectrum for region 1 (top curve, shown in light blue in the electronic edition) includes the effect of pile-up. 
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FIG. 8. — Spatial structure in the brightest regions of G359.23-0.82. The leftmost panel shows the 0.5-8.0 keV image of G359.23-0.82 in the immediate 
vicinity of the peak of X-ray emission; individual 0f'49 X 0f'49 ACIS pixels can be seen. The transfer function in this image is linear, rangin g betw een and 350 
counts per ACIS pixel. The center left panel shows a fit to these data, using a model consisting of two gaussians plus a constant offset (see i|3.4l and Table l2lfor 
details); the transfer function and greyscale range are the same as in the leftmost panel. The center right panel shows the residual between the data and the model, 
again using the same transfer function and greyscale range. The rightmost panel shows the model again, but this time resampled using 0" 12 X 0" 12 pixels, and 
using a linear greyscale ranging between 0.2 and 27.3 counts per 0" 12 X 0" 12 pixel. In the first and f ourth panels, the ellipse shows our revised 3-tr confidence 
limits on the position of PSR J1747-2958, as determined by the radio timing results presented in j|3.5l and in Table|3l 
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-0.1 0.1 0.2 

X (pc) 

FIG. 9. — Axially symmetric hydrodynamic simulation of a bow shock, generated by a pulsar with an isotropic wind moving from right to left through a 
homogeneous ambient medium. The position of the p ulsar is indicated by the "+" symbol; regions referred to in the text are indicated. Simulation parameters 
for the pulsar and ambient medium are as described in SI4.3I The greyscale is on a logarithmic scale; the units in the scale bar are log 10 p, where p is the ambient 
mass density in units of 10~ 24 g cm -3 . 



